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Abstract 

We study the Dirichlet problem in Lipschitz domains and with boundary data in Besov spaces, 
for divergence form strongly elliptic systems of arbitrary order, with bounded, complex-valued 
coefficients. Our main result gives a sharp condition on the local mean oscillations of the coefficients 
of the differential operator and the unit normal to the boundary (which is automatically satisfied 
if these functions belong to space VMO) guaranteeing that the solution operator associated with 
this problem is an isomorphism. 



1 Introduction 

1.1 Formulation of the main result 

A fundamental theme in the theory of partial differential equations, which has profound and intriguing 
connections with many other subareas of analysis, is the well-posedness of various classes of boundary 
value problems under minimal smoothness assumptions on the boundary of the domain and on the 
coefficients of the corresponding differential operator. The main result of this paper is the solution 
of the Dirichlet problem for higher order, strongly elliptic systems in divergence form, with complex- 
valued, bounded, measurable coefficients in Lipschitz domains, and for boundary data in Besov spaces, 
under sharp smoothness assumptions. In order to be more specific we need to introduce some notation. 
Let m, I € N be two fixed integers and, for a bounded Lipschitz domain Q in M. n , n > 2, with 



outward unit normal v = {y\, ...,v n ) consider the Dirichlet problem for the operator 

C(X,D X )U:= Y, D a {A aP {X)D^U), 

\a\=\P\=m 



1.1) 



i.e., 



D a (A a/3 (X)D^U) = 

\a\ = \(3\=m 



for AGO, 



;i.2) 



d k U 



gk on d£l, < k < m — 1. 



(—id/dxi) 



{—id/dx n ) an if a = (ai, a n ). The coefficients A a p are 



Here and elsewhere, D L 

I x I matrix-valued functions with measurable, complex-valued entries, for which there exists some 



*2000 Math Subject Classification. Primary: 35G15, 42B20, 35J55, 35J40. Secondary 35J67, 35E05, 46E39. 
Key words: higher order elliptic systems, Besov spaces, weighted Sobolev spaces, mean oscillations, BMO, VMO, Lipschitz 
domains, Dirichlet problem, traces, extension operators 

The work of authors was supported in part by from NSF, DMS and FRG grants as well as from the Swedish National 
Science Research Council. 



1 



finite constant k > (referred to in the sequel as the ellipticity constant of C) such that 



£ ||4tf 11^(0) < K' 1 (1.3) 
|a| = |/3|=m 

and such that the coercivity condition 

K X E (Atf (*)£>" VW, ^ Q V(X)) > « Yl \\ D ° V HL(n) (1-4) 

a l = l/^l =m |a|=m 

holds for all C'-valued functions V G Co°(f2). Throughout the paper, SRz denotes the real part of 
z G C and (•, •} stands for the canonical inner product in C'. Since, generally speaking, v is merely 
bounded and measurable, care should be exercised when defining iterated normal derivatives. For the 
setting we have in mind it is natural to take d k /dv k := (Y^j=i Cj^/dxj) \^ =v or, more precisely, 

—i^ Y —,v a Ti[D a U], 0<k<m-l, (1.5) 

dv Oil L J> V J 

\a\=k 

where Tr is the boundary trace operator and v a := v" 1 ■ ■ ■ v^ n if a = («i, a n ). Now, if p G (1, oo), 
a G (— 1/p, 1 — 1/p) are fixed and p(X) := dist(X,9f2), a solution for (jl .2f) is sought in W p ' a (Q), 
defined as the space of vector- valued functions for which 

J \D a U(X)\ p p{X) pa dX) lP < oo. (1.6) 

|a|<m 

In particular, as explained later on, the traces in (|1.5p exist in the Besov space Bp(dfl), where s := 
1 — a — 1/p G (0, 1), for any IA G W™' (O). Recall that, with da denoting the area element on dil, 

f€B;(dn)^\\f\\ BSp{dQ) :=\\f\\ Lp{dn) + [[ [ ^'l^Jl da x da Y ) 1/P < oo. (1.7) 

The above definition takes advantage of the Lipschitz manifold structure of d£l. On such manifolds, 
smoothness spaces of index s G (0, 1) can be defined in an intrinsic, invariant fashion by lifting their 
Euclidean counterparts onto the manifold itself via local charts. Nonetheless, the very nature of the 
problem investigated in this paper requires the consideration of higher order smoothness spaces on 
dQ, in which case the above approach is no longer effective. An alternative point of view has been 
introduced by H. Whitney in |63| where he considered higher order Lipschitz spaces on arbitrary 
closed sets (see also C. Fefferman's article [18) for related issues). An extension of this circle of ideas 
pertaining to the full scale of Besov and Sobolev spaces on irregular subsets of M. n can be found in the 
book |30j by A. Jonsson and H. Wallin. Here we further refine this theory in the context of Lipschitz 
domains. For the purpose of this introduction we note that one possible description of these higher 
order Besov spaces on the boundary of a Lipschitz domain 0, C M ra , and for m G N, p G (l,oo), 
s G (0, 1), reads 

B™- 1+S (dn) = the closure of [(i)^D a V|an)| a |<m-i = V G C °°(IR n )} in fl|(02) (1.8) 

(making no notational distinction between a Banach space X and X N = X © • • • X). A formal 
definition, which involves higher order Taylor remainders in place of f(X) — f(Y) in (|1.7p . along with 
other equivalent characterizations of B™~ 1+S (d£l) can be found in §7.1. Given (|1.5p - (|1.6|) and (|1.8|) . a 
necessary condition for the boundary data {gk}o<k<m~i in (|1.2p is that 

9k = 1/01 f°» < fc < m — 1, for some / = {/ a }| Q |< w -l 6 B^~ 1+S {d^). (1.9) 

|a|=Jfc 
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Let BMO and VMO stand, respectively, for the John-Nirenberg space of functions of bounded 
mean oscillations and the Sarason space of functions of vanishing mean oscillations (considered either 
on ft, or on 9ft). Our main result, pertaining to the well-posedness of the problem (|1.2p . then reads 
as follows. 

Theorem 1.1 Let ft C W 1 be a bounded Lipschitz domain whose Lipschitz constant is < M , and 
assume that the operator C in satisfies (jl.3|) - (jl.4p . Then there exists c > 0, depending only on 
M and the ellipticity constant of C, with the following significance. 

Given p G (l,oo) and s G (0,1) seta := 1 — s — l/p. Then the Dirichlet problem (|1.2|) with boundary 
data as in (|1.9j) has a unique solution IA for which (|1.6)) holds provided the coefficient matrices A a p 
and the exterior normal vector v to <9ft satisfy 

dist(z^,VMO(dft)) + dist{A aP ,YMO(n)) < cs(l-s)[p 2 (p-iy 1 + s- 1 (l-s)- 1 ^ .(1.10) 

\a\ = \P\=m 

In particular, the above Dirichlet problem is well-posed for each p G (1, oo) and s G (0, 1) granted that 
A al3 6 VMO(ft) and v G VMO(dft). 

The above result (along with its inhomogeneous version, presented in §8.1) is sharp. See §8.2 for a 
discussion. Note that MBMO(an) +£|a|=|0|=m[A*/3]BMO(«) < cs(l - s)(p 2 (p - + s- l (l - s)~ l y^ 
suffices for (jl.lOp to hold. Other corollaries of independent interest are presented below. 



1.2 Some consequences of the main result and of its proof 

In the proof of Theorem 11.11 we shall actually work with a related condition in place of (|1.10|) , which 
quantifies the local oscillations of the coefficient matrices and the unit normal. To state this formally, 
we need one final piece of terminology. By the BMO mod VMO character of a function F G L\(Q) we 
shall understand the quantity 



{F}^ := limsup ( sup / / F(X) - F(Y) 
e^o \(B s ) n JB £ nn JB s nn 



dXdY 



V 



(1.11) 



where (B £ )n stands for the collection of arbitrary balls centered at points of £1 and of radius e, and 
the barred integral is the mean value. In a similar fashion, if (B E )gQ is the collection of n-dimensional 
balls with centers on d£l and of radius e, and if / G L\{dQ), we set 



{f}*,dn ■= limsup sup 

e->0 \(B £ ) an 



f f 

JB e ndn Jb 



B E ndn 



f(X) - f(Y) 



da x day 



(1.12) 



For an arbitrary function F we obviously have (with the dependence on the domain dropped) {-F}* < 
2dist (F, VMO), where the distance is taken in BMO. Moreover, as a consequence of a result due to 
D. Sarason (cf. Lemma 2 on p. 393 of [50]), there exists C > such that dist(F,VMO) < C{F}*. 
Thus, all together, {F}* ~ dist (F, VMO) so that condition (ll.lOp becomes equivalent to 



Yl i A al3}*,n < cs(l-s)(p 2 (p-l)- 1 + s- 1 (l-s] 

\a\=\/3\=m 



-1 



(1.13) 



A corollary of our main result is that, under the hypotheses of Theorem 1 the problem 

Y D a {A aP {X)D^U) = in ft, 
\a\=\p\=m 



(1.14) 



t z' 7 ' Tr [D^ IA\ = / 7 on dtl, 



I7I < m — 1, 
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has a unique solution hi which satisfies (|1 .6j) whenever 



/ := {/ 7 }| 7 |<m-l G i? p m " 1+S (^)- (1-15) 

However, an advantage of the classical formulation (jl,2p over ()1 . 14j) is that, in the former case, the 
number of boundary conditions is minimal. For a domain Si C K 2 with boundary of class C 1+r , 
i < r < 1, and for real, constant coefficient, scalar operators, the limiting case p = oo of the Dirichlet 
problem (|1.14p has been considered by S. Agmon in |2j. Exploiting the special nature of the layer 
potentials associated with the equation in the two-dimensional setting, he has proved that there exists 
a unique solution hi £ C m_1+S (0), < s < r, whenever / 7 = i^D 7 V\dn, M < m — 1, for some 
function V G C' m ~ 1+S (Q). See also [3 J for a related problem. 

The innovation that allows us to consider boundary data in Besov spaces as in (|1.15p is the 
systematic use of weighted Sobolev spaces such as those associated with the norm in (]1.6p . In relation 
to the standard Besov scale in M n , we would like to point out that, thanks to Theorem 4.1 in [29] on 
the one hand, and Theorem 1.4.2.4 and Theorem 1.4.4.4 in |22j on the other, we have 



a 



1 - a - I E (0, 1 - 1/p) W^ a (n) i? p m - 1+s+1/p (0) 



p 



(1.16) 

a = 1 - s - i G (-l/p,0) => B^ 1+S+1/P (fi) W p m ' a (ft). 

Of course, VF™' a (0) is just the classical Sobolev space W™{£1) when a = 0. 

Remarkably, the classical trace theory for ordinary Sobolev spaces in domains with smooth bound- 
aries turns out to have a most satisfactory analogue in this weighted context and for Lipschitz domains. 
One of our main results in this regard is identifying the correct class of boundary data for higher order 
Dirichlet problems for functions in W p ' a (Q). In the process, we establish that 

U € Wp m ' a (ft) and Pi = on dQ for < k < m - 1 
U belongs to the closure of Cfi°(Q) in W p m,a (Q), 

which provides an answer to the question raised by J.Necas in Problem 4.1 on page 91 of his 1967 
book [47] , In the context of unweighted Sobolev spaces and for smoother domains, such a result has 
been known for a long time (cf., e.g., P. Grisvard, S.M. Nikol'skii and H. Triebel's monographs |22j, 
|48j, |58j and the references therein). 

As a consequence of the trace theory developed in §7 we have that, in the context of Theorem II .1\ 



E 

\a\<m- 



\\Tr[D<*U)\\ B s im) ~ (j2 J p(X)X 1 - 8 )- 1 \D Q U(X)\'>dx\ , (1.18) 
-l \M<™ ^ / 



uniformly in hi satisfying C(X,Dx)hi = in £1. The estimate (I1.18P can be viewed as a significant 
generalization of a well-known characterization of the membership of a function to a Besov space 
in W 1 - 1 in terms of weighted Sobolev norm estimates for its harmonic extension to Wt (see, e.g., 
Proposition T on p. 151 of E. Stein's book [57] ). 



1.3 A brief overview of related work 

Broadly speaking, there are two types of questions pertaining to the well-posedness of the Dirichlet 
problem in a Lipschitz domain f2 for a divergence form, strongly elliptic system (jl.ip of order 2m with 
boundary data in Besov spaces indexed by s and p. 

Question I. Granted that the coefficients of C exhibit a certain amount of smoothness, identify the 
indices p, s for which this boundary value problem is well-posed. 
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Question II. Alternatively, having fixed the indices s and p, characterize the smoothness of d£l and of 
the coefficients of C for which the aforementioned problem is well-posed. 

These, as well as other related issues, have been a driving force behind many exciting, recent devel- 
opments in partial differential equations and allied fields. An authoritative account of their impact is 
given by C.Kenig in [31] where he describes the state of the art in this field of research up to mid 
1990's. One generic problem which falls under the scope of Question I is to determine the optimal 
scale of spaces on which the Dirichlet problem for a strongly elliptic system of order 2m is solvable in 
an arbitrary Lipschitz domain O in W 1 . The most basic case, that of the constant coefficient Laplacian 
in arbitrary Lipschitz domains in IR n , is now well-understood thanks to the work of B. Dahlberg and 
C. Kenig [14J, in the case of L p -data, and D. Jerison and C. Kenig [29], in the case of Besov data. The 
case of (I8.15D for boundary data exhibiting higher regularity (i.e., s > 1) has been recently dealt with 
by V. Maz'ya and T. Shaposhnikova in [41] where optimal smoothness conditions for d£l are found in 
terms of the properties of v as a Sobolev space multiplier. Generalizations of (|8.15p to the case of 
variable-coefficient, single, second order elliptic equations have been obtained in |44] . |45| . [46] . 

In spite of substantial progress in recent years, there remain many basic open questions, particularly 
for I > 1 and/or m > 1, even in the case of constant coefficient operators in Lipschitz domains. In this 
context, one significant problem is to determine the sharp range of p's for which the Dirichlet problem 
for strongly elliptic systems with L p -boundary data is well-posed. In |49j, J.Pipher and G.Verchota 
have developed a L p -theory for real, constant coefficient, higher order systems L = J2\ a \=2mAaD a 
when p is near 2, i.e., 2 — e < p < 2 + e with e > depending on the Lipschitz character of f2. 
On p. 2 of [49] the authors ask whether the L p -Dirichlet problem for these operators is solvable in 
a given Lipschitz domain for p £ (2 — e, + e )> an d a positive answer has been recently given 

by Z. Shen in |51| . Let us also mention here the work [1] of V. Adolfsson and J.Pipher who have 
dealt with the Dirichlet problem for the biharmonic operator in arbitrary Lipschitz domains and with 
data in Besov spaces, |61| where G. Verchota formulates and solves a Neumann-type problem for the 
bi-Laplacian in Lipschitz domains and with boundary data in L2, [43] where the authors treat the 
Dirichlet problem for strongly elliptic systems of second order in an arbitrary Lipschitz subdomain 
Q of a (smooth) Riemannian manifold and with boundary data in B^dQ), when 2 — e < p < 2 + e 
and < s < 1, as well as the paper [32] by V. Kozlov and V. Maz'ya, which contains an explicit 
description of the asymptotic behavior of null-solutions of constant coefficient, higher order, elliptic 
operators near points on the boundary of a domain with a sufficiently small Lipschitz constant. 

A successful strategy for dealing with Question II consists of formulating and solving the analogue 
of the original problem in a standard case, typically when £1 = and C has constant coefficients, 
and then deviating from this most standard setting by allowing perturbations of a certain magnitude. 
A paradigm result in this regard, going back to the work of S.Agmon, A. Doughs, L.Nirenberg and 
V.A. Solonnikov in the 50's and 60's, is that the Dirichlet problem is solvable in the context of Sobolev- 
Besov spaces if dO, is sufficiently smooth and if C has continuous coefficients. The latter requirement 
is an artifact of the method of proof (based on Korn's trick of freezing the coefficients) which requires 
measuring the size of the oscillations of the coefficients in a pointwise sense (as opposed to integral 
sense, as in (jl.lip ). For a version of Question II, corresponding to boundary data selected from 
Iir^o 1 B™~ l - k+s {dQ), optimal results have been obtained by V. Maz'ya and T. Shaposhnikova in [39] . 
In this context, the natural language for describing the smoothness of the domain Q, is that of Sobolev 
space multipliers. 

In the smooth context, problems such as (|1.2|> have been investigated by many authors, including 
L.Garding [ST], M.I. Visik [B2], F.E.Browder [7], S.Agmon, A. Doughs and L.Nirenberg [2], [3J, @], 
V. Solonnikov |55j . [56], L.Hormander [27], G. Grubb and N.J.Kokholm [23] . A related result is 
as follows. If C R n is an arbitrary bounded open set, and g € W™^), m € N, is given, then the 
problem 

A m u = in Q, uelf 2 m (!l), D a (u - g)\ dn = for \a\ < m - 1, (1.19) 
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where the boundary traces are taken in a generalized sense, has a unique solution. Building on some 
earlier work of K. Friedrichs, S.L. Sobolev has considered this problem in |53j in the case when d£l 
consists of a finite union of submanifolds of 1" of arbitrary codimension. This result also appears 
in Sobolev's 1950 monograph [53]. For arbitrary domains, the well-posedness of (|1.19p has been 
established by L.I. Hedberg in |24|, |25j. The issue of continuity of the variational solution for higher- 
order equations at boundary points has been studied by V. Maz'ya in |38j . 

1.4 Comments on the proof of Theorem 11.11 and the layout of the paper 

While the study of boundary value problems for elliptic differential operators with rough coefficients 
goes a long way back (it suffices to point to the connections with Hilbert's 19-th problem and De 
Giorgi-Nash-Moser theory), a lot of attention has been devoted lately to the class of operators with 
coefficients in VMO. Part of the impetus for the recent surge of interest in this particular line of 
work stems from a key observation made by F. Chiarenza, M. Frasca and P. Longo in the early 1990's. 
More specifically, while investigating interior estimates for the solution of a scalar, second-order elliptic 
differential equation of the form CIA = F, these authors have noticed in that U can be related to F 
via a potential theoretic representation formula in which the residual terms are commutators between 
operators of Mikhlin-C alder on- Zygmund type, on the one hand, and operators of multiplication by 
the coefficients of C, on the other hand. This made it possible to control these terms by invoking the 
commutator estimate of Coifman-Rochberg- Weiss (|12j). An alternative method, based on maximal 
operators and good-A inequalities, has been developed by L. Caffarelli and I. Peral in |10j . whereas 
when ft = M. n , an approach based on estimates for the Riesz transforms has been devised by T. Iwaniec 
and C. Sbordone in |28j. Further related results can be found in |5j, [9|, jllj . |17| . |35j. 

Compared to the aforementioned works, our approach is more akin to that of F. Chiarenza and 
collaborators [11 J though there are fundamental differences between solving boundary problems for 
higher order and for second order operators. One difficulty inherently linked with the case m > 1 
arises from the way the norm in (II. 6p behaves under a change of variables x : £1 = {(X',X n ) : X n > 
ip(X')} — > R™ designed to flatten the Lipschitz surface <9f2. When m = 1, a simple bi-Lipschitz changes 
of variables such as f2 3 (X',X n ) ^ (X',X n — ip(X')) 6 will do, but matters are considerable 
more subtle in the case m > 1. In this latter situation, we employ a special global flattening map first 
introduced by J.Necas (in a different context; cf. p. 188 in [47J) and then independently rediscovered 
and/or further adapted to new settings by several authors, including V. Maz'ya and T. Shaposhnikova 
in [39], B. Dahlberg, C. Kenig J. Pipher, E. Stein and G. Verchota (cf. [13J and the discussion in [15J), 
and S. Hofmann and J.Lewis in [26]. Our main novel contribution in this regard is adapting this 
circle of ideas to the context when one seeks pointwise estimates for higher order derivatives of h and 
A := x -1 in terms of [Vc/jJbmcKR™- 1 )- 

Another ingredient of independent interest is deriving estimates for D^DyG(x, y) where G is the 
Green function associated with a constant (complex) coefficient system L{D) of order 2m in the upper 
half space, which are sufficiently well-suited for deriving commutator estimates in the spirit of [12J. 
The methods employed in earlier works are largely based on explicit representation formulas for G(x, y) 
and, hence, cannot be adapted easily to the case of general, non-symmetric, complex coefficient, higher 
order systems. By way of contrast, our approach consists of proving directly that the residual part 
R(x,y) := G(x,y) — — y), where $ is a fundamental solution for L(D), has the property that 
D"DyR(x,y) is a Hardy-type kernel whenever \a\ = |/?| = m. See also [5] for a discussion of the 
difficulties encountered when estimating the residual part R(x, y) in the case when C is not necessarily 
symmetric and has complex coefficients. 

The layout of the paper is as follows. Section 2 contains estimates for the Green function in the 
upper-half space. Section 3 deals with integral operators (of Mikhlin-Calderon-Zygmund and Hardy 
type) as well as commutator estimates on weighted Lebesgue spaces. In the last part of this section we 
also revisit Gagliardo's extension operator and establish estimates in the context of BMO. Section 4 
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contains a discussion of the Dirichlet problem for higher order, variable coefficient, strongly elliptic 
systems in the upper-half space. The adjustments necessary to treat the case of an unbounded domain 
lying above the graph of a Lipschitz function are presented in Section 5, whereas in Section 6 we explain 
how to handle the case of a bounded Lipschitz domain. In Section 7 we study traces and extension 
operators for higher order smoothness spaces on Lipschitz domains. Finally, in Section 8, we deal with 
the inhomogeneous version of f)l -2|) : cf. Theorem 18.11 from which Theorem 11.11 follows . 



2 Green's matrix estimates in the half-space 

2.1 Statement of the main result 

Fix two nonnegative integers m, I and let L(D X ) be a matrix- valued differential operator 

L(D X )= £ A « D *> 

\a\=2m 

where the A a 's are constant I x I matrices with complex entries. Throughout the paper, := 
j- \a\QaiQot2 . . . Qan jf a — ( ai} a2) o; n ) £ N™. Here and elsewhere, N stands for the collection of all 
positive integers and No := N U {0}. Assume that L is strongly elliptic, i.e., there exists k > such 

that H\ a \=ra Pallc^ < K^ 1 and 

M(mV,r))ci>K\Z\ 2m \\v\\h, V£eK n , VrjEC 1 . (2.2) 

In what follows, in order to simplify notations, we shall denote the norms in different finite-dimensional 
real Euclidean spaces by | • | irrespective of their dimensions. Also, quite frequently, we shall make 
no notational distinction between a space of scalar functions, call it X, and the space of vector-valued 
functions (of a fixed, finite dimension) whose components are in X. We denote by F(x) a fundamental 
matrix of the operator L(D X ), i.e., an I x I matrix solution of the system 

L{D x )F{x) = 5{x)k in R n , (2.3) 

where 1\ is the I x I identity matrix and 5 is the Dirac function. We consider the Dirichlet problem 

f L{D x )u = f mm, 



Tr [d j u/dxi] = fj j = 0, 1, . . . ,m - 1, 



on 



pn— I 



(2.4) 



where WL := {x = (x',x n ) : x' £ M ra_1 , x n > 0} and Tr is the boundary trace operator. Hereafter, we 
shall identify cM" with R n_1 in a canonical fashion. 

For each y' G R n_1 we introduce the Poisson matrices Poj • • • > Pm-i for the problem (|2.4p . i.e., the 
solutions of the boundary- value problems 

f L(D x )P J (x,y') = 0I l mWl, 

, Qk \ (2.5) 
{ Pj) ( (*', 0), y' ) = 6 jk S(x' - y')Ii for x' G M. n , < k < m - 1, 

where 5jk is the usual Kronecker symbol and < j < m — 1. The matrix-valued function Pj(x,0') is 
positive homogeneous of degree j + 1 — n, i.e., 

P j (x,tf) = \x\ : > +1 - n P j (x/\x\,<Y), x & M n , (2.6) 

where 0' denotes the origin of W 1 ' 1 . The restriction of Pj{; 0') to the upper half-sphere 5+ 1 is smooth 
and vanishes on the equator along with all of its derivatives up to order m — 1 (see for example, §10.3 
in 1341). Hence, 



x m 



Pj(x,(f)\\ clxl <C u xeR n + , (2.7) 
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and, consequently, 

'm 

\\ p ^^^ c \ x _y^\n^-r * eM +' y'^ n -\ (2.8) 

By G(x, y) we shall denote the Green's matrix of the problem (|2.4p . i.e., the unique solution of the 
boundary- value problem 

f L(D x )G(x,y) = S(x-y)I l for x G M. n , 

( & \ , ( 2 -9) 

(^— Gj ((x',0),y) =0// former- 1 , < j < m - 1, 

where y G is regarded as a parameter. We now introduce the matrix 

R(x, y) := F(x - y) - G(x, y), x, y G R+, (2.10) 
so that, for each fixed y G R™ , 

' L(D X ) R(x, y) = for a; G R™, 

( &> ' J\ (2- 11 ) 

Our goal is to prove the following result. 

Theorem 2.1 For all multi-indices a, (3 of length m 

\\D°DPR(x,y)\\ cM <C\x-y\- n , (2.12) 

for x,y G R™, where y := (y', — y n ) is the reflection of the point y G WL with respect to dW}. 

In the proof of Theorem 12.11 we distinguish two cases, n > 2m and n < 2m, which we shall treat 
separately. Our argument pertaining to the situation when n > 2m is based on the following useful 
estimate for a parameter dependent integral. 

Lemma 2.2 Let a and b be two non-negative numbers and assume that ( G M> N . Then for every e > 
and < 5 < N there exists a constant c(N, e,5) > such that 

dv < c(N,e,5) 

RN (\n\ + a) N +e(\r] -C| + b) N ~ s ~ a £ (|C| + a + b) N ~* ' 1 ' 

The proof is postponed for §2.4, as to prevent disrupting the flow of the presentation. 



-r] ((x',0),y)= (^jF) ((x',0)-y) for x' G R™' -1 , < j < m - 1. 



2.2 Proof of Theorem 12.11 for n > 2m 

In the case when n > 2m there exists a unique fundamental matrix F(x) for the operator (|2.ip which 
is positive homogeneous of degree 2m — n. We shall use the integral representation formula 

R(x,y) = R (x,y) + ... + R m - 1 (x,y), x,yeRl, (2.14) 

where R(x,y) has been introduced in (|2.10p and, with Pj as in (|2.5p . we set 

Rj(x,y):= f Pj(x,?) (®F)((£,0)-y)de, < j < m - 1. (2.15) 

JR"- 1 \OXn / 

Then, thanks to (|2.7p we have 

f x m d£' 

\RAx,y)\vi*i<C \ — " - — -T--TT7T-, — ^ — ( 2 -16) 
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Next, using Lemma 12.21 with 



N = n-l , a = x n , 5 = 2m - j - 1, (2.17) 
e = m-j , b = y n , ( = y' -x', 



we obtain from (|2. 16j) 



ll^-folOllc'*' <7n /i^^^i ^-2m+,- » 0<j<m-l. (2.18) 

(|y' - x'| + x„ + y n ) n 2m+ i 

Summing up over j = 0, . . . , m — 1 gives, by virtue of (I2.14|) . the estimate 

\\R(x,y)\\ cM <C\x-y\ 2m ~ n , x,yGR n + . (2.19) 

To obtain pointwise estimates for derivatives of R(x, y), we make use of the following local estimate 
for a solution of problem (|2.4f) with / = 0. Recall that W£ stands for the classical L p -based Sobolev 
space of order s. The particle loc is used to brand the local versions of these (and other) spaces. 

Lemma 2.3 [see |4j] Let £ and Co be functions inCo°(W l ) such thatQo = 1 in a neighborhood of suppQ. 
Then the solution u G W 2 m (R^, loc) of problem fl23]) with f = and fj G Wp +1 ~ j ~ 1/p {W l ~ 1 , loc), where 
k>m and p G (l,oo), belongs to W~ +1 (M™ , loc) and satisfies the estimate 

m—l 

IIHIh^ +1 (r») - c (Y1 IICo/^II^M-i-j-i/K^n-ij + IICo«llLpCwp)> ( 2 - 2 °) 

where C is independent of u and fj . 

Let B(x,r) denote the ball of radius r > centered at x. 

Corollary 2.4 Assume that u G W™(M?l.,loc) is a solution of problem (|2.4|) with f = and /j G 
Crfc+i-jQfcfi-l j oc ) r /j en f or an y ze wf and p>0, 

m— 1 j 

sup |V fc w|<cfp-* sup M + V J2 p s +i~ k sup |V' S /A (2.21) 

R^nB(z,p) ^ W*nB(z,2p) _q s=0 R^nBfoZp) ' 

where V' s is the gradient of order s in W 1 " 1 . Here C is a constant independent of p, z, u and fj. 

Proof. Given the dilation invariant nature of the estimate we seek, it suffices to assume that p = 1. 
Given <p G C k+1 ~ 3 '(M™" 1 ) supported in I"" 1 n B(z, 2), we observe that, for a suitable 9 G (0, 1), 

k+X-j 

WH^-i-V^) < Cll^Hi ^-^II^H^^ < C Y, sup |V>|. (2.22) 
Also, if p > n, 

sup |Vfc«| < C \\v\\ w k+i (Kn , v (2.23) 

by Sobolev's inequality. Now, (JlTSTj) follows by combining (|2~2"2"|) . HfZ^ with Lemma EH □ 

Given x, y G R+, set p := \x — y\/5 and pick z G cM™ such that |x — z\ = p/2. It follows that for any 
w G K™ PI -B(z, 2p) we have \x — y\ < \x — z\ + \z — w\ + \w — y\ < p/2 + 2p+\w — y\ < \x — y\/2+\w — y\. 
Consequently, \x — y\/2 < \w — y\ for every w G E" PI B(z,2p), so that, ultimately, 

p»- k sup HV^Cw-^llc.xi < ] _^_ 2m+fc , (2.24) 

i«GM n - 1 nS(z,2p) l x ~~ 2/1 
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for each v G No- Granted (|2,19p and our choice of p, we altogether obtain that 

\\D%R(x,y)\\ 0xl <C k \x-y\ 2m ~ n -\ x,yeR n + , Va G N£, |a| < 



(2.25) 



In the following two formulas, it will be convenient to use the notation Rc for the matrix R associated 
with the operator C(D X ) as in (I2.10|) . By Green's formula 



Rc(y,x) = Rc*{x,y) 



x, y G 



(2.26) 



where the superscript star indicates adjunction. In order to estimate mixed partial derivatives, we 
observe that (|2.26|) entails 



(D^R c )(x,y)= (D%R c *)(y,: 



(2.27) 



and remark that C* has properties similar to C. This, in concert with (|2.25|) and the fact that 
\x — y\ = \x — y\ for x,y E WL, yields 



\\D^ y R(x,y)\\ cM ^Cplx-yl 2 ™-^. 
Let us also point out that by formally differentiating (|2.11|) with respect to y we obtain 
' C{D X ) [D^R c {x, y)\ =0 for x 6 E n , 



(2.28) 



^jD^R)((x',0),y) 
dxl 



dxl 



(-DfF) ((x',0)-y), x'e 



nn— 1 



< j < m — 1. 



(2.29) 



With (|2.28p and (12.290 in place of (12. 191) and (I2.11|) . respectively, we can now run the same program 
as above and obtain the estimate 



\D2DgR(x,y)\\ c i*i <C al3 \x-y\ 2m - n -^-W, Va,/?GN . 



(2.30) 



2.3 Proof of Theorem 12.11 for n < 2m 



When n < 2m we shall use the method of descent. To get started, fix an integer N such that N > 2m 
and let (x, z) i— > Q(x,y, z — C) denote the Green matrix with singularity at (y, C) G K n x IR^ - ™ of the 
Dirichlet problem for the operator C{D X ) + (— A z ) m in the iV-dimensional half-space 



:= {(x,z) : z G R N ~ n , x = (x',x n ), x' G I™" 1 , x n > 0}. 
Also, recall that G(x, y) stands for the Green matrix of the problem (|2.4p . 
Lemma 2.5 For all multi-indices a and (3 of order m and for all x and y in 

D«DlG{x', y)= [ D a x D?g{x, y, -Q d(. 

Proof. The strategy is to show that 



(2.31) 



D5D$G(x,y)f p (y)dy 



D2D$g(x,y,-C)dCMy)dy 



(2.32) 



(2.33) 



for each fp G C^°(]R™), from which (|2.32p clearly follows. To justify (I2.33P for a fixed, arbitrary 
fp G Cq°(M™ ), we let u be the unique vector-valued function satisfying D a u G L 2 (M") for all a with 
led = m, and such that 



' C{D x )u = D P x f p inWl 



<9 J u 
dxi 



(x',0) = 



on R n ~\ < j < m - 1. 



(2.34) 
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It is well-known that for each 7 E Nq 

\D^u(x)\ < C 7 \x\ m - n -^ for |x| > 1. (2.35) 



This follows, for instance, from Theorem 6.1.4 |33j combined with Theorem 10.3.2 [34] . Also, as a 
consequence of Green's formula, the solution of the problem (|2.34|) satisfies 

D a x u{x)= [ D«(-D y fG(x,y)fp(y)dy. (2.36) 



We shall now derive yet another integral representation formula for D x u in terms of (derivatives 
of) Q which is similar in spirit to (12. 36ft . Since N > 2m, (|2.30p implies 

\\D%DPg(x,y,-0\\ clxl <c(\x-y\ + \(\r N . (2.37) 

Let us now fix x € R™ , p > and introduce a cut-off function H £ C°°(R N ~ n ) which satisfies H(z) = 1 
for |z| < 1 and H(z) = for \z\ > 2. We may then write 



u{x)= G(x,y,-0 H(C/p)D^f (y) + (-A <: ) m (H(C/p)u(y)) dyd(, (2.38) 
which further implies 

D a x u{x) - [ D«{-D y f g(x, y, -Q H((/p) f p (y) dy d( 



I \\B%D2g(x,y,-0\\ clx i\u(y)D'l(H((/p))\dC. (2.39) 



< c 

|7|=m 

By ()2.35j) and (|2.37p . the expression in the right-hand side of (|2.39p does not exceed 



c P - m [ d( ! dx-yi + icir^iyr-"^ 

J p<\C\<2p JR"- 1 



'p<|C|<2p 

Jf—n—m 



< cp 1 



[ (\y\+pr N \y\ m - n dy = cp- n . 
This estimate, in concert with (12.370 . allows us to obtain, after making p — » 00, that 

D«u(x) =[[ D^(-D y fg(x, y, -() d( fp(y) dy. (2.40) 

JR™ JR N ~ n 

Now (l2~33l follows readily from this and (pOHjh □ 

Having disposed of Lemma 12.51 we are ready to discuss the 

End of Proof of Theorem 12.11 Assume that 2m > n and let N be again an integer such that 
N > 2m. Denote by J-(x,z) the fundamental solution of the operator C(D X ) + (—A z ) m , which is 
positive homogeneous of degree 2m — N and is sing ular at (0,0) G R n x R N ~ n . Then the identity 

D^F{x)=f D^T(x,-()dC (2.41) 

can be established as in the proof of Lemma 12.51 Combining (|2.41|) with Lemma 12.51 we arrive at 

D2D$R(x',y)= [ D2Dln(x,y,-0<K, (2.42) 

where TZ(x, y, z) := g(x, y, z) — T{x — y,z). Consequently, 

\\D^D^n(x,y,-C)\\ cM <C(\x-y\ + \C\)- N . (2.43) 
by (|2.25|) with k = and N in place of n. This estimate, together with (|2.42|) . then yields (|2.12|) . □ 
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2.4 Proof of Lemma EO 

Write J = J\ + J2 where J stands for the integral in the left side of (|2.13j) . whereas J\ and J 2 denote 

the integrals obtained by splitting the domain of integration in J into the ball B a = {r] € : |ry| < a} 
and M. n \B a , respectively. If \(\ < 2a, then 

L w^^ ^LmS?*- (2 ' 44) 

He nc c 

f ca~ N - £ a N lb N ~ & if a<b, 

Ji < < (2.45) 

\ca- N -^ iia>b, K ' 

so that, in particular, 

|£| < 2a => Ji < ca- £ (|C| + a + &)* _JV . (2.46) 
Let us now assume that \Q\ > 2a. Then 

Jl - L wtw^- mr^ - ca ~' m + ' + 6)S " {2A7) 

which is of the right order. As for J72, we write 

j2 - L\ Ba \r,\»+°{\v-C\ + b)»-' = j2 > 1 + J ** (2 - 48) 

where J 2 1, J 22 are obtained by splitting the domain of integration in the above integral into the set 
{rj : \rj\ > max{a,2|£|}} and its complement in M. n \B a . We have 

q < [ ^7 f drj 



|„|> max {a,6,2]C|} \v\ N+£ (\v\ + b) N - 5 A>|„|>max{aA2|C|} \v\ N+£ (\v\ + b) N ~ S 



|r?|>max{a,&,2|C|} |^| 2jY+£ 5 b N 5 J6> |r?| >max{a,6,2| C| } I 7 ?! JV+£ / 

c c 



(a + o+ ICI)^ 6 - 5 a=(a + 5+ ICI)^ 5 
c 

a e (a + 6~+ (CI) 
There remains to estimate the integral 



< — I i I (2-49) 



^2,2 _ / u)V +f n„ ^/-l , u\N-S - ^2,2 + ^2,2 > ( 2 - 50 ) 



\fl„ W JV+£ (I^-CI + &) JV - 5 



where an d ? are obtained by splitting the domain of integration in J 2 .2 into B^/ 2 \B a and its 
complement (relative to B 2 \^\\B a ). On the one hand, 

r (1) < c f dr] < c ( 2 51) 

J ^ - (Id + 6)^-* J Bl(lrABa \ V \»+* -aF(\C\ + a + b)*s ■ [Z ' bL) 

(2) 

On the other hand, whenever \Q > a/2, the integral J 22 1 ' which extends over all 77's such that |r/| > a, 
2|d > [77I > ICI/2, can be estimated as 

7 (2) c f d V c f d£ 



- - icr +e y B2|cl \ Ba (\v-c\+ b) N - s - icr +£ y B4|cl m + 

\ |£|<6 |£|>b / 
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Consequently, 



(2 ) mm{\C\,b} N 
^2,2 - c |^|jv+ e& AT-5 " ^.03) 



Using |£| > a/2 and the obvious inequality 

min{|C|, b} N • max{|C|, ft}" -5 < |C|* o^ 5 (2.54) 



we arrive at 

< ca- £ (\C\ +a + b) s ~ N . (2.55) 

The estimate (I2.55p . along with (|2.5ip and (|2.50[> . gives the upper bound ca~ e (\(\ + a + b) s ~ N for ^2,2- 
Combining this with (|2.49p we obtain the same majorant for Ji which, together with a similar result 
for J\ already obtained, leads to (|2,13|) , The proof of the lemma is therefore complete. 



3 Properties of integral operators in a half-space 

In §3.1 and §3.2 we prove estimates for commutators (and certain commutator-like operators) between 
integral operators in R™ and multiplication operators with functions of bounded mean oscillations, in 
weighted Lebesgue spaces on W). Subsection 3.3 contains BMO and pointwise estimates for extension 
operators from M n_1 onto Wt- Throughout, given two Banach spaces E,F, we let £(E,F) stand for 
the space of bounded linear operators from E into F, and abbreviate £(E) := £(E,E). Also, given 
p G [1, 00], an open set O C M. n and a measurable nonnegative function w on O, we let L p (0, w{x) dx) 
denote the Lebesgue space of (classes of) functions which are p-ih power integrable with respect to 
the weighted measure w(x)dx on O. Finally, following a well-established custom, A{r) ~ B{r) will 
mean that each quantity is bounded by a fixed multiple of the other, uniformly in the parameter r. 



3.1 Kernels with singularities along dW\_ 

Recall that L p (IR™, x^f dx) stands for the weighted Lebesgue space of p-th power integrable functions 
in K" corresponding to the weight w(x) := xlf, x = (x',x n ) G E". 

Proposition 3.1 Let a G K, 1 < p < 00, and assume that Q is a non-negative measurable function 
on {£ = ((', ( n ) G W 1 ^ 1 xl: ( n > —1}, which also satisfies 

I Q(C',Cn-l)C a - 1/p dC<oo. (3.1) 

Then the operator 

Qf(x):=x- n [ Q (y^)f(y)dy, x = (x', x n ) G R%, (3.2) 

J K n \ X n / 

initially defined on functions f G L P (IR™ ) with compact support in Wl_, can be extended by continuity 
to an operator acting from L p (R™, x^ dx) into itself, with the norm satisfying 

\\Q\WL vi Wl, tfdx)) < [ Q(C,(n-l)Cn a - 1/P d(. (3.3) 
+ JR™ 

Proof. Introducing the new variable £ := {x~ (y' — x'),x~ y n ) G WL, we may write 

\Qf(x)\< f Q((',(n-l)\f(x' + x n £,x n ( n )\dC VxGM™. (3.4) 
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Then, by Minkowski's inequality, 

\\Qf\\L p m,^^) < S Q(C',Cn-i)(/ x a n r\f(x' + x n (',x n ( n )\vdxY /P d( 

= Q(C',Cn-l)C a - 1/P dC)ll/llL p (R;,x-dx), (3-5) 

as desired. □ 
Recall that y := (y', -y n ) if y = (y', y n ) G W 1 ' 1 x K. 

Corollary 3.2 Consider 

f log (M + 2 ) 

:= / I /(y) dy, x = (x', x n ) G M" (3.6) 



\x - y 



Then for each 1 < p < oo and each a G (—1/p, 1 — 1/p) t/ie operator R is bounded from L p (W±, x„ dx) 
into itself. Moreover, there exists c{n), independent of a, p and s := 1 — a — 1/p, for which 



c(n)p 2 c(n) 
(pa + l)(p(l - a) - 1) = s(l-s) 



I on ^I'VF _ /o 7 \ 

l K lljC(L p (R",< p (fa:)) S — : TYTZTi r\ TT ~~ 77i ~- IA 



Proof. The result follows from Proposition 13. II with 

._ log(|Cl + 2) 

e tCj — (|£| 2 + 1 )n/2' 

and from the obvious inequality 2|x — y| 2 >|x — y| 2 + x^. □ 
Let us note here that Corollary 13.21 immediately yields the following. 

Corollary 3.3 Consider 

Kf{x) := I -^4- dy, x G R%. (3.9) 

Jmi f - y| 

T/ien /or eac/i 1 < p < oo and a G (—1/p, 1 — 1/p) the operator K is bounded from L P (IR™ , x^f dx) 
into itself In addition, there exists c(n), independent of a, p and s := 1 — a — 1/p, such that 

\\K\\z (Lp(n ^ dx)) < {pa + l){p{l _ a) _ l) ~ ^ZT7y (3-10) 
Recall that the barred integral stands for the mean- value (taken in the integral sense). 



Lemma 3.4 Assume that 1 < p < oo, a G (— 1/p, 1 — 1/p), and recall the operator K introduced in 
(|3.9|) . Further, consider a non-negative, measurable function w defined on R™ and fix a family of balls 
T which form a Whitney covering o/M™. Then the norm of wK as an operator from L p (IR™, Xn ' dx) 
into itself is equivalent to 

sup -fw(y) p dy. (3.11) 



Furthermore, there exists c(n), independent of w, p, a and s := 1 — a — 1/p, such that 

c(n) / f \ Vp 



14 



Proof. Fix / > and denote by \B\ the Euclidean volume of B. Sobolev's embedding theorem allows 
us to write 



\\ K f\\l x{ B) < c N i^r 1 E \ B r /n \\VjKf\\i p{By 

3=0 



VB G .F. 



(3.13) 



Hence, 



/ \x a n w(x)(Kf)(x)\Pdx < c(n) sup -fw(yfdy [ <* V x^VjKffdx. (3.14) 

JR™ BeF JB J~B™ 0<j<l 

Observing that Xn\V j K f\ < c(n)Kf and referring to Corollary 13.31 we arrive at the required upper 



estimate for the norm of wK. The lower estimate is obvious. 



□ 



We momentarily pause in order to collect some definitions and set up basic notation pertaining to 
functions with bounded mean oscillations. Let / be a locally integrable function defined on M. n and 
define the seminorm 



[/]bmo(R") : = sup f f(x) - f f(y) 
B Jb Jb 



dy 



dx, 



(3.15) 



where the supremum is taken over all balls B in W 1 . If / is a locally integrable function on W\, set 



[/]bmo(R«) : = sup + f{x) - f f(y) dy 
(B) JSnffi" JSnR" 



dx, 



(3.16) 



where, this time, the supremum is taken over the collection (B) of all balls B with centers in R". 
Then the following inequalities are straightforward 

[/]bmo(R") < sup f f f(x) - f(y) dxdy < 2 [/] B MO(K™)- (3-17) 

(B) JBnRI JBnR" 



We also record here the equivalence relation 



.) ~ [Ext/] B MO(R")i 



(3.18) 



where Ext / is the extension of / onto M n as an even function in x n . Finally, by BMO(M" ) we denote 
the collection of equivalence classes, mod constants, of functions / on for which [/]bmo(R™ ) < °°- 



Proposition 3.5 Let b € BMO(M™ ) and consider the operator 

\b(x)-b(y)\ 



Tf(x) :-- 



f - y\ 



f(y)dy, xel" + . 



Then for each p G (1, oo) and a G ( — 1/p, 1 — 1/p) 

T : L p (M™, x^dx) — » L p (Rl, x%> dx) 
is a well-defined, bounded operator, such that if s := 1 — a — 1/p then 

c(n) 



Proof. Given x G W]_ and r > 0, we shall use the abbreviations 



v 



(3.19) 



(3.20) 



(3.21) 



b r (x) := f b(y)dy, 
/B(z,r)nRi 



D r (x) := \b(x) -b r (x)\, 



(3.22) 
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and make use of the integral operator 

Sf(x):= [ f(y) d y, xGM™, (3.23) 

JrI \x-y\ n 

as well as its adjoint S* . For each nonnegative, measurable function / on R™ and each x G MIL , 

r/(z) < sf(x) + s*f(x)+ [ IfalMzjgzglM ^y)^ 

< S/(z) + S*/(z) + c(n) [&]bmo(K ? )*V(z), (3-24) 

where K has been introduced in (|3.9p . Making use of Corollary 13.31 we need to estimate only the 
norm of S. Obviously, 

f \t>x„ (x) — b\r-ri\ (x)\ 

Sf(x) < D Xn (x)Kf(x) + / f(y) dy. (3.25) 

Jr™ \x — y\ 

Setting r = \x — y\ and p = x n in the standard inequality 

- b r (x)\ < c(n) log (~ + l) [b] B MO(R™)' (3-26) 

where r > p, we arrive at 

S7(ar) < D Xn (x)Kf(x) + c(n) [6Wo(R«) i?/(z), (3.27) 

where i? is defined in (13. 6p . Let J 7 be a Whitney covering of with open balls. For an arbitrary 
B G J 7 , denote by 5 the radius of -B. By Lemma 13.41 with w{x) := D Xn (x), the norm of the operator 
D Xn (x)K does not exceed 

supf-/ \D Xn (x)\ p dx) < c(n) sup (i\b(x) -b s (x) \ p dx) 1 + c(n) [&] BMO (r™ ) 

£ ^ J B -B G >/ B 

< c(n) [6]bmo(R«), (3-28) 

by the John-Nirenberg inequality. Here we have also used the triangle inequality and the estimate 
()3.26p in order to replace b Xn (x) in the definition of D Xn (x) by b$ (x). The intervening logarithmic 
factor is bounded independently of x since x n is comparable with 5, uniformly for x G B. With this 
estimate in hand, a reference to Corollary 13.21 gives that 

S : LJWL, Xn p dx) -► LJRl, x7 dx) boundedly 

(3.29) 

for each p G (1, oo) and each a G (— 1/p, 1 — 1/p). 

The corresponding estimate for the norm of S is implicit in the above argument. By duality, it follows 
that S* enjoys the same property and, hence, the operator T is bounded on L p (Wl, x^f dx) for each 
p G (l,oo) and a G (— 1/p, 1 — 1/p), thanks to (|3.24|) and Corollary 13.31 The fact that the operator 
norm of T can be estimated in the desired fashion is implicit in the above reasoning. □ 

3.2 Preliminary estimates for singular integrals on weighted Lebesgue spaces 

We need the analogue of Proposition 13.51 for the class of Mikhlin-Calderon-Zygmund singular integral 
operators. Let S 1 " - " 1 stand for the unit sphere in R n and recall that 

Sf(x)=p.v.[ k(x,x-y)f(y)dy, x G R n , (3.30) 

(where p. v. indicates that the integral is taken in the principal value sense) is called a Mikhlin-Calderon- 
Zygmund operator provided the function k : R ra x (M. n \ {0}) — ► M satisfies: 
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(i) k(x, •) G C°°(IR n \ {0}) and, for almost each x G R n , 

max \\D^k(x,z)\\ Loc{R n xS n-i) < oo. (3.31) 

a|<2n 

(ii) k(x, Xz) = \~ n k(x, z) for each z£l™ and each A G R, A > 0; 

(iii) J S n-i k(x,uj) doo = 0, where doj indicates integration with respect to uj G S n . 

It is well-known that the Mikhlin-C alder on- Zygmund operator S and its commutator [S, b] with the 
operator of multiplication by a function b G BMO(M™) are bounded operators in L P (M") for each 
1 < p < oo. Then 

ll<S||£(L p (R™)) < c{n)pp, || [S,b] \\2(L P (RD) < c{n)pp [6] B mo(RJ). (3-32) 

where 1/p + 1/p' = 1 and c(n) depends only on n and the quantity in (|3.3ip . The first estimate in 
()3.32p goes back to the work of A. Calderon and A. Zygmund (see also the comment on p. 22 of |57J 
regarding the dependence on the parameter p of the constants involved). The second estimate in (|3.32p 
was originally proved for convolution type operators by R. Coifman, R. Rochberg and G. Weiss in |12] 
and a standard expansion in spherical harmonics allows to extend this result to the case of operators 
with variable-kernels of the type considered above. 

We wish to extend (|3.32|) to the case when the Lebesgue measure is replaced by Xn P dx, with 
1 < p < oo and a G (— 1/p, 1 — 1/p). Incidentally, a G (— 1/p, 1 — 1/p) corresponds precisely to the 
range of a's for which w(x) := x^f is a weight in Muckenhoupt's A p class, although here we prefer to 
give a direct, elementary proof. 

Proposition 3.6 Retain the above conventions and hypotheses. Then the operator S and its commu- 
tator [S,b] with a function b G BMO(M") are bounded when acting from L P (R™, x^ dx) into itself for 
each p G (1, oo) and a G (—1/p, 1 — 1/p). Then, with s := 1 — a — 1/p and 1/p + 1/p' = 1, 

\\<S\\z(L p (R^,x^dx)) < c{n)(pp + jfYzrpj )' ( 3 - 33 ) 
II [<S,b] \\z( Lr>(R n « Pdx)) < c(n)(pp + 1 ) [&]bmO(b»V (3-34) 



Z(L p (Rl,x a n p dx)) < c{n)[pp + s ^_^ ) [6]bmO(h»)- 

Proof. Let Xj De the characteristic function of the layer 2 J / 2 < x n < 2 1+J / 2 , j = 0, ±1, . . ., so that 
Yljcz Xj = 2. We then write S as the sum S± + S2, where 

Si := - E XjSxk- (3.35) 
b'-fc|<3 

The following chain of inequalities is evident 

\\ S U\\L p {Rl,x a n "dx) < (E / x ^ S i E 



|fe-j|<3 



P \ 1/P 



< c(n)(^yj5( 2 X ^ a / 2 /)(x)px) 1/P . (3.36) 
i K + i|<3 

In concert with the first estimate in (|3.32j) . this entails 

< c(n)pp'(j2[ (E X k 2 ja/2 \f\) P dx) 1/P 

j \k-j\<3 

< c(n)pp'(f \f{x)\^x^dx) 1/P , (3.37) 

\ /iron / 
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which is further equivalent to 

\\Si\\z(L p {mi,x a n p dx)) <c(n)pp'. (3.38) 
Applying the same argument to [Si , b] and referring to (|3.32p , we arrive at 

II [ s ii b ] \\z(L p (Ri,x a n p dx)) < c(n)pp' [6] B mo(R™)- (3-39) 

It remains to obtain the analogues of (|3.38|) and (|3.39p with 52 in place of <Si. One can check 
directly that the modulus of the kernel of 52 does not exceed c(n) |x — y\~ n and that the modulus of 
the kernel of [52,6] is majorized by c(n) \b{x) — b{y)\ \x — y\~ n - Then the desired conclusions follow 
from Corollary 13.31 and Proposition 13.51 □ 

3.3 BMO type estimates for Gagliardo's extension operator 

Here we shall revisit a certain operator T, extending functions defined on R n_1 into functions defined 
on R™, first introduced by E. Gagliardo in [20j. Fix a smooth, radial, decreasing, even, non-negative 
function £ in R™ _1 such that = for |i| > 1 and 



J C(t) dt = 1. (3.40) 



(A standard choice is £(t) := cexp (—1/(1 — |t| 2 )+) for a suitable c.) Following |20j we then define 

(Tip){x',x n ) := f C(t)(p(x' + x n t)dt, (x',x n )eWL, (3.41) 



acting on functions tp from LiiW 1 1 ,loc). To get started, we note that 

V x ,(T<p)(x',x n ) = J C(t)V<p(x' + tx n )dt, (3.42) 

(Ty){x',x n ) = / ((t)tX7p{x' + tx n )dt, (3.43) 



dx_ 

and, hence, we have the estimate 

IIVzCT^Hi^QRn) <c\\V x ,(p\\ Loo( jg, n -iy (3.44) 

Refinements of (|3.44|) are contained in the Lemmas I3.7H3.8I below. 

Lemma 3.7 (i) For each multi-indices a with \a\ > 1 there exists c > such that 

D2(T<p)(x)\ < c4-HpvV] B MO(R»-i), Vx = (x',x n ) e Rl. (3.45) 

(ii) There exists c > such that 

(Tip)(x) - (p(x') < ci n [V^BMO{R»-i)) Vx = (x', x n ) G . (3.46) 

Proof. Rewriting (|3.43p as 

-£-(T<p)tf,x n )=x 1 r r n [ C(— ) — (V^(e)-/ V^(z)dz)di (3.47) 

OX n J \ X n / X n \ J\z—x'\<x n ' 
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we obtain 



d 

£) 2^-( T ¥')(^) <c\ l7l N BM o(r-i) 



for every non-zero multi- index 7. Furthermore, for i = 1, . . . n — 1, by (13,421) 



9 

— Vj;'(T(/?)(x) = x 



Hence, once again 



5 



(3.48) 



(3.49) 



(3.50) 



and the estimate claimed in (i) follows. Finally, (ii) is a simple consequence of (i) and the fact that 

(?V)Ir™-i = </>• D 

Remark. In concert with Theorem 2 on p. 62-63 in |57J, formula (|3.42p yields the pointwise estimate 



|V(ZV)(x)| < cM{Vtp){x'), x = (x',x n ) G 



11, 



(3.51) 



where M. is the classical Hardy-Littlewood maximal function (cf., e.g., Chapter I in |57J). As for 
higher order derivatives, an inspection of the above proof reveals that 



D a x {T^)(x) < cxlr^iVtpfix'), (x',x n ) G 



(3.52) 



holds for each multi-index a with \a\ > 1, where is the Fefferman-Stein sharp maximal function 

(cf. m) 



Lemma 3.8 IfV x "p G BMO^™" 1 ) then V(I» G BMO(K™) and 

[V(T^)]bmo(k?) < C^'^BMOjl"- 1 )- 



(3.53) 



Proof. Since (T<p)(x', x n ) is even with respect to x n , it suffices to estimate [V x (7V)] BMO (Kn). Let Q r 
denote a cube with side-length r centered at the point rj = (r/,7/ n ) G M n_1 x R. Also let Q' r be the 
projection of Q r on M n_1 . Clearly, 



V x '(T(p)(x',x n ) - V X '<p(x') = x 



l-n 
n 



c 



C-x' 



(v^O-v^'))^. 



(3.54) 



Suppose that \rj n \ < 2r and write 

/ ^x'{T(p)(x', X n ) - V x iip(x') 

JO.r 



Therefore, for \q n \ < 2r 

+ + \X7 x >Tip(x) -X7 y >T<p(y)\dxdy 

J Qr J Qr 



dx < cr 



2-n 



\V<p(£) ~ Vip(z)\dzd£. 



< C^fV^BMOCR"" 1 )- 



< 



2-1 \V x >T<p(x) - V<p(a/)\dx+ f f \Vip{x')-V<p{y')\dx'dy' <c[Vcp] 

JQr J Q'r J Q'r 



(3.55) 



(3.56) 



BMO(R n_1 )- 
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Next, consider the case when \r) n \ > 2r and let x and y be arbitrary points in Q r (rj). Then, using the 

generic abbreviation Je '■= 4- /, we may write 

Je 



\V x ,T<p(x)-V y ,Tip{y)\ < / x^ n C 



l-n/-(£ x '\ „.l-nf(£ V 



2 1 r; n | 



< 



c r 



vy(0 - v^ Q 



2 1 In 



f S|»jnl 

< C [V^BMOCR™- 1 )- 



Consequently, for \r/ n \ > 2r, 

i i \V x >T(p(x) - V y >T(p(y)\dxdy < c [V^BMOfR"- 1 ) 
which, together with (|3.56p . gives 

This inequality and (|3.48|) with \ j\ = imply (|3.53j) . 

4 The Dirichlet problem in for variable coefficient systems 

4.1 The setup 

For 

1 1 

1 < p < oo, — < a < 1 and m € N, 

p p 

we let V p m ' a (Wl) denote the weighted Sobolev space associated with the norm 



(3.57) 
(3.58) 

(3.59) 
□ 



Mlv p m ' a (]R™) : = 



V f \x^- m D p u{x)\ p x p n a dx 



1/p 



(4.1) 



(4.2) 



\(3\<m 



It is easily proved that C^°(1R™) is dense in V^ m ' a (M"). Moreover, by the one-dimensional Hardy's 
inequality (see, for instance, [37J, formula (1.3/1)), we have (with s := 1 — a — 1/p, as usual), 



\U\\\r m , a 



VjT' a (]R™) 



\P\=m 



The dual of V r 



(Rl) will be denoted by Vp 
Consider now the operator 



r) •x v n a dx) IV for u£C °°(i;; 



), where 1/p + 1/p' = 1. 



(4.3) 



L(x,D x )u:= D^A a p{x)x\^\-^Diu) 
\a\,\P\<m 



(4.4) 



where A a p are C' x '-valued functions in Loo(IR+). We shall use the notation L(x, D x ) for the principal 
part of L(x, D x ), i.e., 

L(x, D x )u := D^{A a p{x)D^ x u). (4.5) 

\a\=\p\=m 
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4.2 Solvability and regularity result 

Lemma 4.1 Assume that there exists n = const > such that the coercivity condition 



9ff 



/ Yl (A af }{x)DPu{x), D a u(x)) cl dx > k V \\D^u\\ 2 L , wl) , (4.6) 

+ |a| = |/3|=m |7|=m 

ZioWs /or a// n G C£°(]R+), and i/iai 

X] P^llioo^) < ( 4 - 7 ) 

|a| = |/3|=m 

(i) Let p G (1, oo), a G (—1/p, 1 — 1/p), and suppose that 

g ( l _ J\ Y \\ A ap\\L^(R^ + Y WBMO(iy < ^ ( 4 -8) 

^ ' |a| + ||8|<2m | a | = |fl| =m 

0<|a|,|/3|<m 

where s := 1 — a — 1/p, 1/p + = 1, and (5 satisfies 

pp' H — r) <5 < c(n, to, k) (4-9) 



S (l- S ) 

tuiia a sufficiently small constant c(n, m, k) > 0. T/ien 



L = L(x, D x ) : V p m ' a (R n + ) — ► V~ m ' a {Rl) isomorphically. (4.10) 



(ii) Let pi G (1, oo) and —1/pi < < 1 — l/p%, where i = 1,2. Suppose that |^.g[ ) aoZds mi/j 
Pi and Sj := 1 — a« — 1/pj in place of p and s. Then, if the function u G ^^' ai (IR.!L) is such that 
Lu G V p - m ' ai (M™ ) n Vp~ m > a2 (H£), it /oMouw iaai « G Vp™' a2 (P n ^ 



Proof. The fact that the operator in (|4.10|) is continuous is obvious. Also, the existence of a bounded 
inverse L^ 1 for p = 2 and a = follows from (|4.6[) and (|4.8|) - (|4.9|) with p = 2, a = 0, which allow 
us to implement the Lax-Milgram lemma. We shall use the notation L y for the operator L(y,D x ), 
corresponding to (|4.5p in which the coefficients have been frozen at y G Wt, and the notation G y for 
the solution operator for the Dirichlet problem for L y in WL with homogeneous boundary conditions. 
Next, given u G V p m ' a {Rl), set / := Lu G Vp~ m ' a (M!}:) so that 



+> 



L(x, D)u = f in 
^V,0)=0 onE"- 1 , 0<i <m-l. 

We may now write 

u{x) = {G y f){x)-{G y {L-L y )u){x)-{Gy{L-L)u){x), x G E", (4.12) 

and aim to use (|4.12[) in order to express u in terms of / (cf. (|4.26|) - (|4.27|) below) via integral operators 
whose norms we can control. First, we claim that whenever py| = m, the norm of the operator 

V p m ' a (Rl) 3«h D2(G y (L - L y )u){x) G L^-tf* dy) (4.13) 

x=y 

does not exceed ^ 

+ (1 _ J E MMbmo(R!P- (4-14) 

^ ' | a |=|/3|=m 
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Given the hypotheses under which we operate, the expression (j4,14p is therefore small if 5 is small. 

In what follows, we denote by G y (x,z) the integral kernel of G y and integrate by parts in order 
to move derivatives of the form D" with \a\ = m from (L — L y )u onto G y (x,z) (the absence of 
boundary terms is due to the fact that G y (x, •) satisfies homogeneous Dirichlet boundary conditions). 
That (I4.14p bounds the norm of (14.130 can now be seen by combining Theorem 12.11 with (I3.2ip and 
Proposition 13.61 Let 7 and a be multi-indices with I7I = m, \a\ < m and consider the assignment 



C °°r + )9^fe/ G y (x,z)D c z 



*(*) 
m— \a\ 



dz 



x=y 



After integrating by parts, the action of this operator can be rewritten in the form 

-1 d \m-\a\ 



X 



dz r , 



(-D z ) a G y (x,z) T a (z)dz 



x=y 



(4.15) 



(4.16) 



where 



T a (z) :-- 



if \a 
(_i)«H a l 



m — a 



IV. 



m, 



(t 



\m-\ 



im-\a\ 



dt, 



if lad < m. 



(4.17) 



Using Theorem 12.11 along with ()3.2ip and Proposition 13.61 we may therefore conclude that 



DZ 



-1 d \ m ~\ 



i dz n 



(-D z ) a G y (x,z) T a {z)dz 



x=y 



L P (WL,yn 



On the other hand, Hardy's inequality gives 



< C 



c 



( / 1 



773^) J 1 1 r ° 1 1 L v (K + , rt p dx) ■ 



a II Lj,Qa$,xl p dx) ^ ^^ L p( K l^ x n P dx) 



(4.18) 



(4.19) 



and, hence, the operator (|4.15p can be extended to a linear mapping from Co°(M") to L P (E™ , x^f dx) 
with norm < \pp' + a (i_ s ) j- Next, given u G V™' a (Wl), we let * = ^ a/3 in (|4TT5|) with 

(4.20) 



**p{z) := zW- m A a pD?u(z), \a\ + \0\ < 2m, 
and conclude that the norm of the operator 

V p m > a (R n + ) 3u^DJ (G y (L - L)u)(x) G L p (Wl,y%>dy) 

does not exceed 



x=y 



c 



1-s 



VP + 



1 



a/3 1 1 Loo (MJ)- 



a| + |/3|<2m 
|a|,|/3|<m 



8(1-8). 

It is well-known (cf. (1.1.10/6) on p. 22 of [37]) that any u G V p m,a (l 

u = K{D a u} W \ =m 
where K is a linear operator with the property that 

D a K : L p (M™ , x a P dx) — ► L p (R™ , < p dx) 



(4.21) 



(4.22) 



can be represented as 

(4.23) 

(4.24) 
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is bounded for every multi- index a with \a\ = m. In particular, by (|4.3p . 

\\K{D a U }\ cr \ =m \\ V ^,a {W u ) < C S~ 1 \\{D a u}\ a \ =m \\ Lp{W n >x ap dx y (4.25) 

At this stage, we transform the identity (|4.12p as follows. First, we express the two u's occurring 
inside the Green operator G y in the left-hand side of (I4.12|) as in (I4.23p . Second, for each 7 G No with 
I7I = m, we apply D 7 to both sides of (|4.12p and, finally, set x = y. The resulting identity reads 

{D-y U }\j\=m + S{D°u} wl=m = Qf, (4.26) 

where Q is a bounded operator from Vp~ m ' a (W^ ) into L P (R™ , Xn dx) and S is a linear operator mapping 
L p (R^, xTdx) into itself. Furthermore, on account of P~T3>P~T3]) . (|4T2"T|) - (j4~2"2"]) and (g^SJ), we can 
bound ||S , || £ ( Lp ( R n ja ,«P d:r)) by 

[4/3]bmo(r^)+ s(1 1 _ s) Yl HA*/?|U~(K»))- (4.27) 

V 7 |a| = |/3|= m V ^ | Q | + |0|<2m 

0<|a|,|/3|<m 

Owing to (|4.8p - (|4.9p and with the integral representation formula (|4.26p and the bound (|4.27p in 
hand, a Neumann series argument and standard functional analysis allow us to simultaneously settle 
the claims (i) and (ii) in the statement of the lemma. □ 

5 The Dirichlet problem in a special Lipschitz domain 

In this section as well as in subsequent ones, we shall work with an unbounded domain of the form 

G = {X = (X\ X n ) G M n : X' € R n ~ x , X n > cp(X')}, (5.1) 
where ip : R™- 1 -> E is a Lipschitz function. 

5.1 The space BMO(G) 

The space of functions of bounded mean oscillations in G can be introduced in a similar fashion to 
the case G = M™ . Specifically, a locally integrable function on G belongs to the space BMO(G) if 

1/]bmo(G) == sup / f(X) - I f(Y) dY dX < 00, (5.2) 

(B) JBC\G JBnG 

where the supremum is taken over all balls B centered at points in G. Much as before, 

[/]bmo(G) ~ sup / / f(X) - f(Y) dXdY. (5.3) 
(S) JBnG JBnG 

This implies the equivalence relation 

[/]bmo(G) ~ [/ ^]bmo(k«) ( 5 - 4 ) 
for each bi-Lipschitz diffeomorphism A of onto G. As consequences of definitions, we also have 



]J /j]BMO(G) < c H/llxLoG) [/]BMO(Gf) J where / = (/l> M> ( 5 ' 5 ) 
[/ _1 ]bMO(G) < cll/^Hl^^^BMOCG)- (5-6) 
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5.2 A bi-Lipschitz map A : — ► G and its inverse 

Let if : R™ -1 — > R be the Lipschitz function whose graph is and set M := || V^H^^n-i). Next, 
let T be the extension operator defined as in (|3,4ip and, for a fixed, sufficiently large constant C > 0, 
consider the Lipschitz mapping 

A : IR" 9 (x',x n ) i — > (X',AC n ) € G (5.7) 

defined by the equalities 

X':=x', I n :=CMx n + (T^)(i',x n ) (5.8) 
(see [39j, §6.5.1 and an earlier, less accessible, reference |40j). The Jacobi matrix of A is given by 

V = ( v^czv) cm + d(T<p)/dx n ) ' (5 - 9) 

where 7 is the identity (n — 1) X (n — l)-matrix. Since \d(Ttp)/dx n \ < cM by (j3.43p . it follows that 
det A' > (C - c) M > 0. Thanks to (ET351) and (^7f |) -(^g ]) we have 

X n - <p(X') ~ x n . (5.10) 

Also, based on (|3.53p we may write 

[A']bmo(r™) < c [VvIbmocr™- 1 ) (5.11) 

and further, by (pTl4"j) and (pTl5j) . 

|| J D Q A'(x)|| R „xn < c(M)x-W[V^]bmooii«-i) s Va : |q| > 1. (5.12) 

Next, by closely mimicking the proof of Proposition 2.6 from |41j it is possible to show the existence 
of the inverse Lipschitz mapping x := A -1 : G — > . Owing to (15. 4p . the inequality (|5.1ip implies 

[A' O X] BM0 ( G ) < C [V(y9] BMO (Rn-l). (5.13) 

Furthermore, (I5,12p is equivalent to 

\\(D a \')(x(X))\\ MnXn < c(M,a)(X n - V 9(X / ))- |a| [V^BMO(R"- 1 ), (5-14) 
whenever \a\ > 0. Since x' = (A' o x) _1 we obtain from (|5.6p and (|5.13p 

MbmO(G) < cfV^BMOfR"" 1 )- ( 5 - 15 ) 

On the other hand, using x' = (A' o x) _1 and (|5.14|) one can prove by induction on the order of 
differentiation that, for all X £ G and a £ No with |a| > 0, 

||D<Vp0|| R »x„ < c{M,a) (X n - ^(X , ))- |q| [Vv3]bmo(r- 1) - (5-16) 

5.3 The space V™' a (G) 

Analogously to V^ m ' a (M™), we define the weighted Sobolev space V™' a (G) naturally associated with 
the norm 

\P\\v^{G) ■■= ( E / G l(^-^(^)) l7hm ^W| P (^n-^')r^) 1/P - (5-17) 
|7|<m 

Replacing the function X n — (p(X') by either p(X) := dist (X, dG), or by the so-called regularized 
distance function p Teg (X) (defined as on pp. 170-171 of |57j), yields equivalent norms on V™' a (G). 
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A standard localization argument involving a cut-off function vanishing near dG (for example, take 
r/(p reg /e) where r/ £ Cq°(R) satisfies rj(t) = for |i| < 1 and rj(t) = 1 for \t\ > 2) shows that 

C °°(G) F p m ' a (G) densely. (5.18) 

Next, we observe that for each U € Cg°(G), 

c*iiwii^»(G)<(E / w)r(^-^^')r^) 1/?, < iiwii^.- (G)> (5.i9) 

\y\=m G 

where, as before, s = 1 — a — 1/p. Indeed, for each multi-index 7 with [7I < m, the one-dimensional 
Hardy's inequality gives 

/ \{X n - V {X')p- m D^U{X)\v (X n - ip{X')f a dX 
Jg 

<{C/s) P Y. I \D a U(X)f(X n -ip(X')r a dX, (5.20) 

a|=m 

and the first inequality in (|5.19)) follows readily from it. Also, the second inequality in (|5.19p is a 
trivial consequence of (|5.17p . Going further, we aim to establish that 

Cl |M|y™. a (]R*-) < \\U O X-|| V m,a (G) < C 2 IMIv"'"^™)) (5-21) 

where c\ and c 2 do not depend on p and s, and x : G — > is the map introduced in §5.2. Clearly, 
it suffices to prove the upper estimate for \\u o x:\\y™' a (Q) in (|5.21|) . As a preliminary matter, we note 
that 

IF {<*&))) = {{>c'*{X)i)l =D u)(x(X)) 

n 

+ Y, (D T u)(x(X))Y,c„HllD^x t (X), (5.22) 

l<|r|<|7| cr j=l j 

where 

er = (try)) E ^' = 7 ' l ^'' - lj E^'l ~ X ) = ItI ~ l r l- ( 5 - 23 ) 
In turn, (|5.22p - (|5.23p and (|5.16p allow us to conclude that 

\D^ (u(x(X)))\ < c x ] n hhl \D T u(x)\, (5.24) 

i<|r|<M 

which, in view of (|5,10p . yields the desired conclusion. Finally, if, as usual, p' = p/(p — 1), we set 

V- m > a {G):=(v r ;>- a (G))\ (5.25) 

5.4 Solvability and regularity result for the Dirichlet problem in the domain G 

Let us consider the differential operator 

CU = C(X,D X )U= D a (^(X)D f3 U), XeG, (5.26) 

\a\ = \P\=m 
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whose matrix-valued coefficients satisfy 

IISMWG) < (5-27) 

|a| = |/3|=m 

This operator generates the sesquilinear form £(•, ■) : ^, m ' a (G) x T^ ) T' _a (G') — > C where p' is the 
conjugate exponent of p, defined by 

C(U,V):= [ {X a p(X)DPU(X), D a V{X))dX. (5.28) 

\ a \ = \f3\= m JG 

We assume that 

RC{U,U)>k \\D^ U\\l 2{G) holds for all U G V™'°(G). (5.29) 

\f\=m 



Lemma 5.1 (i) Let p € (1, oo), 1/p + 1/p' = 1, — 1/p < a < 1 — 1/p, and set s := 1 — a — 1/p. Also, 
suppose that 

[V^bmocm™- 1 ) + [^]bmo(g) < S, (5.30) 

\a\ = \f3\=m 

where 5 satisfies 

{ pp> + s (i _ gj ) s (i _ g) < C '^ n ' m ' K ' II V V 9 IIl 00 (iR"- 1 )) ( 5 - 31 ) 
TOi/i a sufficiently small constant C, independent of p and s. Then 

C(X,D X ) : V™' a {G) — ► V- rn ' a {G) isomorphically. (5.32) 

(ii) Suppose that pi € (l,oo) ; —1/pi < ai < 1 — l/pi, i = 1,2. and that (|4.9I) holds with pi and 
Si := 1 — at — 1/pi in place of p and s. If U € V^ ,ai (G) and CU £ Vp^ m ' ai (G) n Vp 2 m,a2 (G), then 

Proof. We shall extensively use the flattening mapping A and its inverse studied in §5.2. The 
assertions (i) and (ii) will follow directly from Lemma 14.11 as soon as we show that the operator 

L(U o A) := (CU) o X (5.33) 

satisfies all the hypotheses in that lemma. The sesquilinear form corresponding to the operator L will 
be denoted by L(u,v). Set u{x) := U(X(x)), v(x) := V(X(x)) and note that (I5.22|) implies 

D?U(X) = ((*>* (\(x))t;)l =D u)(x)+ Yl K^x^D^uix), (5.34) 

1<M<|/3| 

D a V(X) = ((x'*(X(x))Of =D v)(x)+ J2 K aT (x)x^-^D^v(x), (5.35) 

l<|r|<|a| 

where, thanks to (I5.16|) . the coefficients K~ T satisfy 

II-^ttIUooOk^.) - c [V^bmocr™- 1 )- (5.36) 
Plugging (|5.34p and (|5.35p into the definition of C(U,V), we arrive at 

C(U,V) = L (u,v) + £ [ (A aP {x)x^ + ^- 2m D p u{x), D a v(x))dx, (5.37) 

l<|c*|,|/3|<m 
\a\ + \0\<2m 



26 



where 



L (u,v)= I ((^oX)((K'*o\)Ol =D u,((x / *oX)^ =D v)detX'dx. (5.38) 

\ a \ = \(3\=m Jw + 

It follows from (|5.34p - (|5.36p that the coefficient matrices A a p obey 

^ \\A a ^\\ Loo (Mn_) < ckT 1 [V^] BMO (R"-i), (5.39) 



l<\ct\,\P\<m 
\a\ + \/3\<2m 



where c depends on m, n, and ||V(/3|| ioo ( K n-i). We can write the form Lq(u,v) as 

V / (A aP (x)D l3 u(x), D a v(x))dx, (5.40) 

|a]=|/?l=™ + 

where the coefficient matrices A a p are given by 

A a(3 = detX> Yl ^(^'°A)(2l 7T oA), (5.41) 

|7| = |r|=m 

for some scalar homogeneous polynomials Pj^ of the elements of the matrix x'(A(x)) of degree 2m. 
In view of (l531)- (l5~T5|) . 

^ [^]bmo(rj) < cl^^V^BMOtr- 1 ) + ^ Wbmo(g)). ( 5 -42) 

a| = |/3|=m |ct| = |/3|=m 

where c depends on n, m, and HV^Hi^nRn-i). By (|5.39l) 

\L(u,u) - L (u,u)\ < c5||u|[^ m , 0(M „ ) (5.43) 

and, therefore, 

KL (u,u) > ®£(U,U)-c6\\u\\ 2 vr ,o my (5.44) 
Using (|5.29p and the equivalence (cf. the discussion in §5.3) 

||^|| y m,0 (G) ~ IMIy™.°(]Rn); ( 5 - 45 ) 

we arrive at (14. 6p . Thus, all conditions of Lemma 14.11 hold and the result follows. □ 



6 The Dirichlet problem in a bounded Lipschitz domain 
6.1 Background 

Let f2 be a bounded Lipschitz domain in R n which means (cf. [57] , p. 189) that there exists a finite 
open covering {£^}i<j<Ar of 8ft with the property that, for every j G {1, N}, Oj n O coincides with 
the portion of Oj lying in the over-graph of a Lipschitz function ipj : M n_1 — > R (where R n_1 x M is 
a new system of coordinates obtained from the original one via a rigid motion). We then define the 
Lipschitz constant of a bounded Lipschitz domain Q C M n as 

inf (maxlHV^II^^-i) :l<j<N}), (6.1) 

where the infimum is taken over all possible families {(Pj}i<j<N as above. Such domains are called 
minimally smooth in E. Stein's book |57J. It is a classical result that, for a Lipschitz domain Q, the 
surface measure da is well-defined on dfl and an outward pointing normal vector v exists a.e. on d£l. 
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0<\f3\<m' 



We denote by p{X) the distance from X G R n to dSl and, for p, a and m as in (|4.ip . introduce the 
weighted Sobolev space lp m ' a (r2) naturally associated with the norm 

[ \p{X)W~ m D p U{X)\ p p{Xf a dX) lh . (6.2) 

<\P\<m 

One can check the equivalence of the norms 

\\^\\v^' a (Q) ~ l|Preg W llv™.°(n)' ( 6 - 3 ) 

where /9 reg (^) stands for the regularized distance from X to (in the sense of Theorem 2, p. 171 in 
|57J). Much as with (I5.18p . it is also easily proved that Cq°(0) is dense in V™' a (Q) and that 



V p m ' a (f2) ^ 



J \D^U(X)\ p p(X) pa dXy /P (6.4) 



\P\=m' 

uniformly for U G Cq°(Q). As in (|5.25p . we set 

v- m > a (n)-.= (v p 7>- a (n)y. (6.5) 

Let us fix a Cartesian coordinates system and consider the differential operator 

AU = A(X,D X )U := Yj D a (A a( s(X)D (i U), lefi, (6.6) 

\a\ = \/3\=m 

with measurable I x I matrix-valued coefficients. The corresponding sesquilinear form will be denoted 
by A(U, V). Similarly to (|5.27p and (|5.29p we impose the conditions 

Y H-^opllwn) ^ ( 6 - 7 ) 

\a\=\f)\=m 

and 

^A{U,U) > k Y \\ D1 U\\l 2 (n) for a11 u e IC' ^)- ( 6 - 8 ) 

|7|=m 

6.2 Interior regularity of solutions 

Lemma 6.1 Let O C 1" be a bounded Lipschitz domain. Pick two functions Ti, Z G C^°(r2) such that 
HZ = TL, and assume that, for a sufficiently small constant c(m,n, k) > ; 

Y Wbmo(JI) < ( 6 - 9 ) 

\a\ = \0\=m 

where 

S < (6.10) 
pp' 

IfU G % m (0, loc) for a certain q < p and AU £ W~ m {Q, loc), then U G W p m (ft, Zoc) and 

l|WW|| W m(n) < C (WHA(-,D)U\\ w - m ^ + \\ZU\\ w ^(q)). (6.11) 
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Proof. We shall use the notation Ay for the operator A(Y, Dx), where 7 6!! and the notation <E>y 
for a fundamental solution of Ay in W l . Then, with star denoting the convolution product, 

HU + $y *(A-A Y )(HU) = $y*(HAU) + $y*([A,H}(ZU)) (6.12) 

and, consequently, for each multi-index 7, [ — ^ | = m, 

D^{HU)+ D a+ ^Y*((A a f3-A af 3(Y))DP(HU)) 

\a\ = \(3\=m 

= D 7 $ y * (HAU) + D n/ $ Y * {[A,H]{ZU)). (6.13) 
Writing this equation at the point Y and using (I3.32|) . we obtain 

(l-Cpp'5) W{HU)\\ Lpm 

\~l\=m 

< C( P ,K)(\\HAU\\ w - mm + \\[A,H](ZU)\\ w - m{n) ). (6.14) 
Let p' < n. We have for every V G W™(Q), the closure of C^°(Q) in W™(Q), 

[ ([A,H](ZU),V)dX =\A(HZU,V)-A(ZU,HV)\ 
Jn 

< ciWZUWw^wWVWw^in) + \\ZU\\ W rn n (f^llVH^m-l (6.15) 

^ n-\-p _p_n 

n — p' 

By Sobolev's theorem 

\\ zu \\w™-\n) - c \\ zu \\w^ n (n) (6-16) 

n+p 

and 

\ W m-l (n) < C \\V\\ W rn {n y (6.17) 



p' n 



Therefore, 



/ ([AW](2W),V)dX < c ||Z^|| wm ( n)l|V|| W7( n) (6.18) 
Jn htp p 

which is equivalent to the inequality 

\\[A,H}(ZU)\\ Wp - m{a) <c\\ZU\\ wmpn (n) . (6.19) 



n+p 



In the case p' > n, the same argument leads to a similar inequality, where pn/{n + p) is replaced by 
1 + e with an arbitrary e > for p' > n and e = for p' = n. Now, (16. lip follows from (|6.14p if p' > n 
and p' < n, q > pn/(n + p). In the remaining case the goal is achieved by iterating this argument 
finitely many times. □ 

Corollary 6.2 Let p > 2 and suppose that §6M) and (|6~TUj) hold. If U G VT 2 m (0,Zoc) and AU G 
VF-' m (ft, Zoc), i/ien U G % m (^, Zoc) and 

\\HU\\wp m < C(\\ZA(;D)U\\ w - m{n) + ||ZW|| wr -i (n) ). (6.20) 

Proof. Let Zq denote a real- valued function in C^°(fi) such that TCZq = TL and Z$Z = Zq. By (|6.1ip 

ll^ w llw™(n) < C(ll^-4(-)^ ) )^llw- p - m (^) + \\ 2 oU\\w?, n (n)) (6-21) 
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and it follows from (|6,8p that 

\\Z U\\ 2 wr{n) < cK- l $tA(Z U,Z U). (6.22) 

Furthermore, 

\A{Z U,Z U) -A(U,Z$U)\ KcK^WZUWy^-i^ \\Z U\\ w ™ m . (6.23) 

Hence 

W^o^Ww^in) - CK -AMWwz™^) u Ww^(n) + K_ ll-^llw™- 1 ^) \\ 2 o^\\w^(n)) ( 6 -24) 
and, therefore, 

\\Z U\\ Wi n {n) < cK^iWZ AU\\ w - m{n) +K~ 1 ||^W|| w ,m-i (n) ). (6.25) 

Combining this inequality with (|6.2ip we arrive at (|6.20p . □ 

6.3 Invertibility of A : V™' a (Q) — ► V p - m ' a (Q) 
Recall the quantities defined in (jl.lljl - (jl.!2p . 

Theorem 6.3 Let 1 < p < oo, < s < 1, and set a = 1 — s — 1/p. Furthermore, let Vt be a bounded 
Lipschitz domain in W 1 . Suppose that the differential operator A is as in §6.1 and that, in addition, 

{Atfkn + M*,dn < S, (6.26) 

\a\=\0\=m 

where 

(pp'+ //On ( 6 - 27 ) 
V s{± — s)J s{± — s) 

for a sufficiently small constant c > independent of p and s. Then the operator 

A : V™> a {Q) — ► V- m ' a {Vt) isomorphically. (6.28) 

Proof. We shall proceed in a series of steps starting with 

(i) The construction of the auxiliary domain G and operator C. Let e be small enough so that 

V / / \A a p{X)-A a p(Y)\dXdY< 25 (6.29) 
I, ,„i JB r nnJB r nfi 

\a\=\p\=m 

for all balls in {B r }^ with radii r < e and 

/ / v(X) - v(Y) daxdo-y < 25 (6.30) 

for all balls in {B t }q^ with radii r < e. We fix a ball B £ in (B £ )qq and assume without loss of 
generality that, in a suitable system of Cartesian coordinates, 

nnB £ = {X = (X', X n ) E B £ : X n > <p{X')} (6.31) 

for some Lipschitz function <p : W 1 " 1 — > M. Consider now the unique cube Q(e) (relative to this 
system of coordinates) which is inscribed in B £ and denote its projection onto M n_1 by Q'(e). Since 
V(p = —v'/v n , it follows from (|6.30p that 



JBL JBL 



V(f(X')-V(p(Y') dX'dY' <c(n)5, (6.32) 
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where B' r = B r n W 1 1 , r < e. Let us retain the notation </j for the mirror extension of the function tp 
from Q'(e) onto R n_1 . We extend A a p from Q(e) n O onto Q(e)\£l by setting 

A*/3p0 := A a/3 (X', -X n + 2^(X')), X G Q(e)\fi, (6.33) 

and we shall use the notation 2l Qj g for the periodic extension of A a p from Q(e) onto M n . 

Consistent with the earlier discussion in Section 5, we shall denote the special Lipschitz domain 
{X = (X',X n ) : X' G R n-1 , X n > (p(X')} by G. One can easily see that, owing to the 2en~ 1/2 - 
periodicity of <p and A a p, 

^2 [^Ibmo(G) + [V^bmoCM"- 1 ) < c(n)S. (6.34) 

\a\ = \P\=m 

With the operator A(X, Dx) in Q, we associate the auxiliary operator C(X, Dx) in G given by (|5.26j) . 

(ii) Uniqueness. Assuming that U G V™' a {Q) satisfies LIA = in $7, we shall show that U G 
y™i°(Qy This will imply that IA = which proves the injectivity of the operator in (|6.28j) . To this 
end, pick a function H G Cq > (Q(e)) and write C{TLU) = [£, TL\U. Also, fix a small > and select 
a smooth function A on R;L, which is identically 1 on [0, 1] and which vanishes identically on (2, oo). 
Then by (ii) in Lemma 15-H 

C(HU) - [C, H] (A(p ieg /0)U) G V 2 ~ m '°(G) n y~ m ' a (G). (6.35) 

Note that the operator 

[£, H]p-g = ^ m '°( G ) — V m ' a ( G ) ( 6 - 36 ) 

is bounded and that the norm of the multiplier p reg A(p reg /0) in V™' a (G) is O{0). Moreover, the same 
is true for p = 2 and a = 0. The inclusion (|6.35|) can be written in the form 

C(HU) + M(ZU) G V~ m ' a (G) n F 2 ~ m '°(G), (6.37) 

where Z G C^°(M n ), ZH = H and is a linear operator mapping 

V™' a (G) V p - m ' a (G) and V 2 m '°(G) V 2 mfl (G) (6.38) 

with both norms of order O(0). 

Select a finite covering of 0, by cubes Qj(e) and let {Hj} be a smooth partition of unity subordinate 
to {Qj(s}}. Also, let Zj G C^(Qj(e)) be such that Wj-Zj = Hj. By we denote the special Lipschitz 
domain generated by the cube Qj(e) as in part (i) of the present proof. The corresponding operators 
£ and M will be denoted by Cj and M.j, respectively. It follows from (I6.37P that 

Hj U + Ys^i 1 M i Z i Z k){H k U) G V^ a (fl) n V™'°(n). (6.39) 

k 

Taking into account that the norms of the matrix operator Cj Aij Zj Z^ in the spaces ^^"'"(Q) and 
V™' u (n) are O(0), we may take 9 > small enough and obtain HjU G V^(ft), i.e., U G V 2 m ' u (n). 
Therefore, C : V™' a (Q) — ► VJ, _m ' a (f2) is injective. 

(iii) ^4 priori estimate. Let p > 2 and assume that U G VJ, m,a (f2). Referring to Corollary 16.21 and 
arguing as in part (ii) of the present proof, we arrive at the equation 

Hj U + XX t1 M i Z i Zk ^ Hk U ) = T > ( 6 - 4 °) 

k 

whose right-hand side satisfies 

ll^"llvs»-(n) < c(\\AU\\ v - m ,a {n) + (6.41) 
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for some domain oj with w C !!. Since the V^ m,a (f2)-norm of the sum in (|6.40p does not exceed 
C6\\U\\ym,a^^ we obtain the estimate 

Mv^w < c(\\Au\\ Vp - m>a{Q) + \P\\wr\^- ( 6 - 42 ) 

(iv) End of proof. Let p > 2. The range of the operator A : V™' a ($l) — > V p ~ m ' a {Q) is closed by 
(|6.40p and the compactness of the restriction operator: V p m ' a {VL) — > W^"~ 1 (u;). Since the coefficients 
of the adjoint operator C* satisfy the same conditions as those of £, we may conclude that C* : 
V^' a {^l) — > V^,, m ' a (^) is injective. Therefore, C : V^ m,a (0) — > V p m ' a (J7) is surjective. Being also 
injective, C is isomorphic if p > 2. Thus £* is isomorphic for p' < 2 so £ is isomorphic for p >2. □ 



7 Traces and extensions 

7.1 Higher order Besov spaces on Lipschitz surfaces 

Let C R n be a bounded Lipschitz domain and, for m 6 N, 1 < p < oo and —1/p < a < 1 — 1/p, 
consider a new space, W™' a (£l), consisting of functions U G /oc) with the property that p a D a U £ 

L p {Vt) for all multi-indices a with |a| = m. We equip W^'^O) with the norm 

|o|=m 



where a; is an open non-empty domain, u C f2. An equivalent norm is given by the expression in (jl.6p . 
We omit the standard proof of the fact that 

C°°(Tt) ^ W™< a {VL) densely. (7.2) 

Recall that for p G (1, oo) and s G (0, 1) the Besov space Bp(dQ) is then defined via the requirement 
(|l,7p . The nature of our problem requires that we work with Besov spaces (defined on Lipschitz 
boundaries) which exhibit a higher order of smoothness. In accordance with |30j, |57J, |63j, we now 
make the following definition. 

Definition 7.1 For p G (1, oo), m G N and s G (0,1), define the (higher order) Besov space 
B p n ^ l+S {d^l) as the collection of all families f = {fa}\a\<m-i °f measurable functions defined on 
d£l, such that if 

R a (X,Y) := f a (X) - ^f* + p(Y)(X -Yf, X,Yedn, (7.3) 

\(3\<m-l-\a\ 

for each multi-index a of length < m — 1, then 

:= \\f<x\\L P (dn) (74) 

|«|<m— 1 



+ v (/ / . — 'fr'YTn, .aW* 

yJanJanlX -YI?'- m - 1+ °- , ° B+ "- 1 

\a\<rn—l 



< OO. 



It is standard to prove that B™ +s (di~l) is a Banach space. Also, trivially, for any constant k > 0, 

Y Wfa\\L p (dn)+ Y (// \ X - Yt^ls-nHn-l ^xdcTy) '* (7.5) 

\a\<m-l |a|<m-l x ,redn 

\X-Y\<k 
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is an equivalent norm on B™~ 1+S (dft). 

A few notational conventions which will occasionally simplify the presentation are as follows. Given 
a family of functions / = {/ a }| Q |< m _i on dO, and X E W 1 , Y, Z E dQ, set 

P a (X,Y) := Yl ^f a+ a(Y)(X -Yf, Va:|a|<m-1, (7.6) 
\/3\<m-l-\a\ 

P(X,Y) := P (0j ... j0) (X,y). (7.7) 

Then 

R a (Y,Z)=f a (Y)-P a (Y,Z), Va:\a\<m-1, (7.8) 
and the following elementary identities hold for each multi- index a of length < m — 1: 

iWl%P a (X,Y) = P a+0 (X,Y), \/3\ <m-l-\a\, (7.9) 

P a (X,Y)-P a (X,Z) = J^a+^Z^X-Yf. (7.10) 

\P\<m-l-\a\ 

See, e.g., p. 177 in [57J for the last formula. 

We now discuss how 2?™~ 1+s (<90) behaves under multiplication by a smooth function with compact 
support. 

Lemma 7.1 For each p E (1, oo), m E N and s E (0, 1), the Besov space B™~ 1+S (dft) is a module 
over Cft°(K n ), granted that for each f = {/ Q }| a |< m -i 6 B™~ 1+s {dVL) and ij) E Cg°(M n ) we set 

Proof. Let / = {/4 H < m _i E £™- 1+s (<9f}) and V € Cg°(R n ) be arbitrary and set 

xj} ■= {i "I Tr [D a f[\ . (7.12) 

I J |a|<m— 1 

Denote by R a (X,Y) the remainder (17. 3p written for ip f in place of / = {/a}| a i< TO -i- Also, let 
P a (X, Y") and P(X, Y) be the polynomials defined in (17. 6D . (|7.7f) with the components of / replaced 
by those of ip. 

Next, fix some a E Nq with \a\ < m — 1 and, for each 1,7 6 d£l write 

iUx,y)= E -^""^W/rW 

" E fi[ E { ^T^D^{Y)f T {Y)](X-Yf. (7.13) 

5|<m— 1— 1« <T+r=ct+<5 

The crux of the matter is establishing that for a fixed k > there exists C > independent of / such 
that 

I./? (a" y)p 

\ X - y | P ("m-i+ a -|a|)+n-i ^ ^H/H ^-i+^ny ^ 7 - 14 ) 

|X-r|<« 
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To get started, we note that the first sum in (|7.13p can be further expanded as 

/3+7=a /3+7=a 

+ E 

+ E ljhh(X>Y)Pj(X>Y). (7.15) 

Now, if \X — Y\ < k, the first sum in the right-hand side of (|7.15p is pointwise dominated by 
C\X — y| m ~l Q l £\ i <m _ 1 |/ 7 (X)| hence, when raised to the p-th power and then multiplied by \X — 
y|-p(™-i+s-M)+«-i ; it is dominated by 

C E // < x yj?J p{1 - a) *'x<toY<C E U^W (7-16) 

\-y\<m-l x ,Yedfi h\<m-l 

\X-Y\<k 

which suits our goal. Similarly, for \X — Y\ < k, the second sum in the right-hand side of (|7.15p is 
— C X^|7 |<|q 1^7(^5 ^)l thus, as before, its contribution in the context of estimating the left-hand side 
of (|7.14p does not exceed 

c E // nf Jfi^-HHn-! * c WfW V-w (7 - 17) 

M<H x,Y£dn 

\X-Y\< K 

As for the last sum in the right-hand side of (17.151) . we employ (|7.9p . (17. 7p . Leibniz's rule and the 
definitions of P(X,Y), P(X,Y), in order to successively transform this sum into 

E y- ] i m D P x P(X,Y)i^D^P(X,Y)=^D^[P(X,Y)P(X,Y)} 



iHD x f E -^-^D^{Y)f T {Y){X - Y) a+T I . (7.18) 



<J\T 

|cr|,|r|<m— 1 



Note that D%[(X - Y) a+T ] = unless it is possible to select 5 G such that a + r = a + 5. In 
the latter situation, we use i^D^[(X — Y) a+S ] = ("+^)- — Y) s and re-write the last expression in 
(FTTHD as 



|5|<2(m-l)-|a!| cr+T=a+<5 



^[ E ^P^ kl ^W)/.(E)](^-^) 5 - (7-19) 



Now, the second sum in (|7.13p cancels the portion from (|7.19p corresponding to the case when 
\5\ < m — 1 — \a\, and the remaining terms in this sum are < C\X — Y| m ~l a Yl\j\< m—i |/ 7 (X)|. 
Consequently, in the context of (|7.14p . their contribution is estimated as we did in (|7.16p . 

This analysis establishes (|7.14[). Since, trivially, \\(^f) a \\L p (dQ) < c 12\j\<m-l Wf-r\\L p (dU), the proof 
of the lemma is finished. □ 

Typically, the previous lemma is used to localize functions / G B™~~ 1+S (dfl) in such a way that the 
supports of their components are contained in suitably small open subsets of <9f2, where the boundary 
can be described as a graph of a real- valued Lipschitz function defined in R n_1 . Such an argument, 
involving a smooth partition of unity, is standard and will often be used tacitly hereafter. 

We next discuss a special case of the general trace result we have in mind. 
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Lemma 7.2 For each 1 < p < oo, < a < 1 — 1/p and s := 1 — a — 1/p, the trace operator 



Tr : W^ a (n) — » B*(dn) (7.20) 

zs well-defined, linear, bounded, onto and has Vp' a (Q) as its null-space. Furthermore, there exists a 
linear, continuous mapping 

£:B s p (dn)^W^ a (Q), (7.21) 
called extension operator, such that Tr o £ = I (i.e., a bounded, linear right-inverse of trace). 

Proof. By a standard argument involving a smooth partition of unity it suffices to deal with the case 
when Q is the domain lying above the graph of a Lipschitz function ip : R n_1 — > R. Composing with 
the bi-Lipschitz homeomorphism Wf 3 (X',X n ) i— > (X',cp(X') + X n ) £ f2 further reduces matters 
to the case when Q = Wf, in which situation the claims in the lemma have been proved in (greater 
generality) by S.V. Uspenskh in |59j (a paper preceded by the significant work of E. Gagliardo in |20j 
in the unweighted case). □ 

We need to establish an analogue of Lemma 17.21 for higher smoothness spaces. While for SI = Wf 
this has been done by S.V. Uspenskh in |59j, the flattening argument used in Lemma 17.21 is no longer 
effective in this context. Let us also mention here that a result similar in spirit, valid for any Lipschitz 
domain f2 but with B m ~ l+s+l l p (Vl) in place of W™' a (p) has been proved by A. Jonsson and H. Wallin 
in [30j (in fact, in this latter context, these authors have dealt with much more general sets than 
Lipschitz domains). The result which serves our purposes is as follows. 

Proposition 7.3 For 1 < p < oo, — 1/p < a < 1 — 1/p, s := 1 — a — 1/p £ (0, 1) and m G N, define 
the higher order trace operator 

tr m _! : W^ a (Q) — (7.22) 

by setting 

tr m _i U:= \i\ a \Tr[D a U}\ , (7.23) 

l J \a\<m— 1 

where the traces in the right-hand side are taken in the sense of Lemma 7.2. Then (|7.22j) - (|7.23|) is a 
well-defined, linear, bounded operator, which is onto and has V p m,a (U) as its null-space. Moreover, it 
has a bounded, linear right-inverse, i.e., there exists a linear, continuous operator 

£ : B™- l+s {dSl) — ► W™' a (n) (7.24) 

such that 

f = {f a } H < m ^i^Bp- 1+s (dn)^i\ a \Tr[D a (£f)] = f a , V«: \a\ < m - 1. (7.25) 

In order to facilitate the exposition, we isolate a couple of preliminary results prior to the proof of 
Proposition 17.31 The first is analogous to a Taylor remainder formula proved by H. Whitney in |64j 
using a different set of compatibility conditions than (|7.26p below. 

Lemma 7.4 Assume that ip : M™- 1 -»■ R is a Lipschitz function and define $ : R n 1 — » dfl =— > M. n 
by setting $(X') := (X',ip(X')) at each X' € R n_1 . Next, consider the special Lipschitz domain 
£1 := {X = (X',X n ) € R n : X n > ip(X')} and, for some fixed m G N, a system of sufficiently nice 
functions {/a}| a |<m-i with the property that 

f) n 

[f a ($(X'))] = J2U+eM X ')) d k®j( X % l<fc<n-l, \a\<m-2, (7.26) 



dX k 



35 



where {ej}j is the canonical orthonormal basis in R n . Then, with R a (X,Y) defined as in (|7.3p . the 
following identity holds for each multi-index a of length < m — 2 and a.e. X',Y' € K ri ~ . 

R a ^{x')MY')) = Y, E - { [fa+^mr + tix' -¥>)))- / Q+7+e .($(y')) 

j=l [ 7 [=m-2-|a| 7 ' 

x($(x') -$(y' + t(x' -y'^v^y' + ^x'-y')) ■ (x'-y')dt. ( 7 - 27 ) 

Proof. We shall prove that for any system of functions {fa}\a\<m-l which satisfies (|7.26p . any multi- 
index a € Nq with \a\ < m — 2 and any I £ N with I < m — 1 — |a|, there holds 



/ a ($(x')) - £ l/ a+/3 (*(r'))(*(*') - $(y')/ 



l/3|<' 



E E 



1 



7! Jo L 



/ a+7+e .(d>(y' + t(x' - y'))) - / Q+7+e .($(y')) 



i=i | 7 |=/-i 

x($(x') - $(y' + i(x' - y'))) 7 v$ j (y / + t(x' - y')) • (x' - y') dt. (7.28) 

Clearly, (f73D) follows from j73]) and (I7T28]) by taking 2 := m - 1 - \a\. 

In order to justify (17.281) we proceed by induction on I. Concretely, when I = 1 we may write, 
based on (|7.26p and the Fundamental Theorem of Calculus, 



UHx')) - f a (HY')) - ^/ Q+ej (d>(y'))(^(x') - *,(y')) 



(7.29) 



1 d_ 



M^Y' + t(X' - Y>)))}dt-J2fa + eM Y ')) i j t \®j{Y' + t{X'-Y')) 



3=1 



dt 



n , rl 



= E{ / [/«+e,W + - y'))) - / a+e .($(y'))] v^(y' + *(x' - y')) • (x' - y') dt}, 

as wanted. 

To prove the version of (|7.28|) when I is replaced by / + 1 we split the sum in the left-hand side 
of (|7.28p . written for I + 1 in place of I, according to whether \j3\ < I or |/3| = I + 1 and denote 
the expressions created in this fashion by S\ and S% respectively. Next, based on (|7.26p and the 
Fundamental Theorem of Calculus, we write 



/ a+7+e .($(y' + t(X' - Y'))) - f a+1+e ^{Y ! )) 
n ft 



(7.30) 



/ f a+1+ej+ e k (HY' + T (x> - y')))v$ fe (y' + r(x' - y')) • (x> - y>) dr 
k =i Jo 



and use the induction hypothesis to conclude that 

-l r t 



s i=E E I?/ / We J+ e fc ($(y' + r(x'-y')))(<i>(X)-^(y , + t(X-y'))r 

j,fc=l | 7 |=Z-1 ~' ^° ^° 

xV$j(7' + t{X' - Y')) ■ (X 1 - y')V$ fc (y' + r(X' - Y')) ■ (X' - Y') dr dt. (7.31) 
Thus, if we set 

F j {t):=* j {X')-* j {Y' + t{X'-Y% l<J<n, (7.32) 
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and use an elementary identity to the effect that for any M n -valued function F = (F\, ...,F n ), 

1 ^ 1 



d_ 

dt 



(7.33) 



we may express S\ in the form 



n fi ft 

si = -EE / [fa +P+ e k mY' + t(x> - y'))) - f a+p+ek mr))} 

k=im=l Jo Jo 



dt 



F(tf V$ fc (E' + t{X' - Y')) ■ (X' - Y') drdt 



k=l \fi\=l 



1 rt 



JO 



dt 



±F(tf]FUr)dTdt. 



(7.34) 



Note that after changing the order of integration and using the Fundamental Theorem of Calculus, 
the first double sum above corresponds precisely to the expression in the right-hand side of (|7.28|) 
written for I + 1 in place of I. By once again changing the order of integration and relying on (|7.33p . it 
becomes apparent that the second double sum in (|7.34|) is — S^- Thus, Si + 52 matches the right-hand 
side of ([725$ with I replaced by I + 1, proving (j7T28l) . □ 

Corollary 7.5 Under the assumptions of Lemma \7.4[ for each multi-index a of length < m — 2 the 
following estimate holds 

|7|=m— 1 

where the constant C depends only on n, p, s and \\^^p\\L 00 (R n - 1 )- 
Proof. The identity (fLZTJ) gives 

|i? a ($(X'),$(y'))l 

< (7|x' — y'l™- 1 -!"! 



(7.36) 



£ f 1 |/ 7 (*(y' + t(x' - F'))) - / 7 (*(r' 

_ 1 -/o 



))!<« 



|7|=m— 1 

for each X',Y' G M n_1 , where the constant C depends only on n and HV'i&ll^nKn-i) which, in turn, 
is controlled in terms of || VyHifRn-i). If we now integrate the p-th power of both sides in (|7.36|) 
for X' ,Y' G R n_1 , use Fubini's Theorem and make the change of variables Z' := Y' + t(X' — Y'), we 
obtain 



\R a (X,Y)\i 



-da x da Y <c! [ \R a mX>)MY W 



on Jan \X - y|j»(tn-i-M+.)+ 



|+s)+n-l 



dX'dY' 



^ E 

|7|=m— 1 

^ E 

| 7 |=m— 1 



«-l ./IBn-l 



^\f^{Y' + t{X'-Y')))-f^{Y'W 



«-l ./IBi-l 



_ y/lps+n— 1 



W _ \ps+n— 1 



eft dAW 



< 



° E 

| 7 |=m— 1 



'7iiB|(an)' 



(7.37) 
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since |V$(X')| ~ 1 and \§(X') - $(Y')\ ~ |X' - Y'\, uniformly for X',Y' G W 1 - 1 . □ 

After this preamble, we are in a position to present the 

Proof of Proposition 17.31 We divide the proof into a series of steps, starting with 

Step I: The well-definiteness of trace. Thanks to (|7.2p . it suffices to study the action of the trace 
operator of a function U G C°°(f2) W™' a (£l). If we now set 

f a :=i lal Tr[D a U], Va:|a|<m-1, (7.38) 

it follows from Lemma 17.21 that these trace functions are well-defined and, in fact, 

^ ll/a||s|(an) < C\\U\\ w ^,a^ n y (7.39) 

\ct\<m— 1 

In order to prove that / := {/ a }|a|<m-i belongs to B™~ l+s (dQ), let R a (X,Y) be as in (|7.3jl . Our 
goal is to show that for every multi- index a with \a\ < m — 1, 

(f I i ^tt^tti ^ i ^ x ^yV /P < C , l|Z/|| w m.. fnv (7.40) 

To this end, we first observe that if |a| = m — 1 then the expression in the left-hand side of (j7.40H 
is majorized by c(j m J m \f a {X) - f a (Y)\P\X - y|-d»+n-i) daxdayY P which, by fTB^ . is indeed 

< C|| To treat the case when |a| < m — 2 we assume that is locally represented 
as {X : X n > (p(X')} for some Lipschitz function tp : M n_1 — > M and, as before, set $(-X~') := 
(X',(p(X')), X G M™ -1 . Then (17361) holds, thanks to (17381) . for every multi-index a of length 

< m — 2. Consequently, Corollary 17.51 applies and, in concert with (|7.39p . yields (I7,40p , This proves 
that the operator (|7.22[) - ()7.23f) is well-defined and bounded. 

Step II: The extension operator. We introduce a co-boundary operator 6 which acts on / = 
{/a}|o|<m~i G Bp- l+s {dty according to 

(£f)(X)= [ K.(X,Y)P(X,Y)da Y , (7.41) 
Jan 

where P(X,Y) is the polynomial associated with / as in (|7.7H . The integral kernel /C is assumed to 
satisfy 

/ K{X, Y) day = 1 for all (7.42) 

\D X K,{X,Y)\ <c aP {X) 1 ~ n -\° l \ VXen,VYedn, (7.43) 
where a is an arbitrary multi-index, and 

K(X,Y) = Q if \X-Y\ > 2p{X). (7.44) 
One can take, for instance, the kernel 

x-y \ f r ( x- z N 



JC(X,Y):= V [ ^ r, i — da z , (7.45) 

where r] G C^°(i?2), ?7 = 1 on B±, rj > and x is a positive constant depending on the Lipschitz 
constant of dQ. Here, as before, p Teg (X) stands for the regularized distance from X to <9f2. 
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For each X € ft, and Z € dO, and for every multi- index 7 with \ j\ = m we then obtain 

unsf{x)= £ i / DWx,r)(P^x,r)-P£(x,z))da y . (7.46) 

|c|>l 

Fix ^ > 1 and denote by B(X,R) the ball of radius R centered at X. We may then estimate 

\unsf{X)\ p < c V pixy^-f \p (x,y) - p p (x, z)\p day 

J B(x,fi P (x))ndn 



c*+/3=7 
|a|>l 



^ C H E P(Xr pH -f \R 5+ p(Y,Z)\ p \X-Y\^da Y , 
a+P=n w|+|a|< m _i 7 B(x, w (x))nai 

|a|>l 

< C V p{X) p ^ rn ^-f \R T (Y,Z)\ p da Y , (7.47) 

| T |< m _l B(X,w(X))ndSl 

where we have used Holder's inequality and (|7.10p . Averaging the extreme terms in (|7.47p for Z in 
B(X, fip(X)) PI dfl, we arrive at 

\D~<£f(X)\ p <C p(X) p( ^- m ^ 2(n -^ j f \R T {Y,Z)\ p da Y da z . (7.48) 



[t|<to-1 Y.zeen 

\X-Y\,\X-Z\<iip(X) 

For each multi-index 7 of length m we may then estimate: 

I iD^SfiX^piX)^ 1 -^- 1 dX 
Jn 

^ C J2 I I \ R r(Xi Z )\ P ( [ p{X) pi -- m+1 - s+ \ T ^- 2n+1 dX) d<J Y d 

1 r. -, Jan Jan ^ J ' 

\T\<m-l X6 n 

\X-Y\,\X-Z\<np(X) 

\T\<m— 1 

<C\\f\f B ^ +S{my (7.49) 
by (|7.4p . where we have used the readily verified fact that there exists C > such that 

J p ( X y{-m+ls+\T\)-2n+l dx < jnrjy _ Z | P (- m+ l+| r |- s )_„+l } 

x en 

|x-y|,|x-z|< M p(^) 



\X-Y\,\X-Z\<^p(X) 

for any <9fl and r £ Nq with |r| < m — 1. This proves that the operator (|7.4ip is well-defined 

and bounded in the context of (|7.24|) . 

Step III: The right-invertibility property. We shall now show that the operator (I7.4ip is a right-inverse 
for the trace operator (|7.22|) . i.e., whenever / = {/ 7 }| 7 |< m _i £ B™~ 1+s (dQ), there holds 

/ 7 = »N Tr[L»^/] (7.51) 
for every multi-index 7 of length < m — 1. To this end, for I7I < m — 1 we write 

WSf(X)-£ y f(X)= i2 -W\ I D x lC(X,Y)(Pp(XX)-Pp(X,Z))da Y , (7.52) 



H>1 
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where 



fiy/PQ := / /C(X,F)P 7 (X,F)d(ry, lea (7.53) 
Jan 

Estimating the right-hand side in (|7.52p in the same way as we did with the right-hand side of (|7.46p . 
we obtain using the boundedness of dfl 

|r|<m— 1 



< C ll/H ^- 1+S(an) - (7.54) 



In a similar fashion, we check that 



/ \v(wsf(x) - sjixwpixy-v*- 1 dx 

Jan 



<C Y I I ; ^ T iyy,^ 1 dcydcz < C Wff+m-i+s^- (7.55) 



|r|<m— 1 



an Jan 



The two last inequalities imply D~ t £f — £f £ Vp ,a (d) and, therefore, 

Tr {D~<£f - E 1 f) = 0. (7.56) 

Going further, let us set 

Eg(X):=f JC{X, Y) g(Y) day, X € fi. (7.57) 
Jan 

A simpler version of the reasoning in Step II yields that E maps Bp(dfl) boundedly into Wp' a (Q). 
Also, a standard argument based on the Poisson kernel-like behavior of )C(X, Y) shows that Tr Eg = g 
for each g £ B*(dSl). Based on (175)1 -(177 )1 and (1733)) we have 

\£ 7 f(X) - Ef y (X)\P + p(Xy\V(£j(X) - Ef^XW 

<C P( X ) m { \f 1+ p{Y)\Vday. (7.58) 

|/3|< m -i-M J B(X,fip(X))nda 

I/3|>1 

Consequently, for an arbitrary Whitney cube Q C of side- length I, we have 

/ \£,f(X)-Ef,(X)\Pp(X)-^- 1 dX+ [ \V{£ 1 f{X)-Ef 1 {XWp{Xf-^- l dX 

JQ JQ 

< C l pm - s) [ \f, +l 3(Y)\ P d<TY, (7.59) 

l /3 |< m _i_l 7 l JaanxQ 
\p\>i 

where xQ is the concentric dilate of Q by some fixed factor x > 1. Summing over all cubes Q of a 
Whitney decomposition of 0, we find 

\\£~ff - Ef^Wyi.a^ < C ^ H/aJUptsn) (7.60) 

|ct|<m— 1 

which implies 

Tr(£j-Ef 7 ) = 0. (7.61) 
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Finally, combining (|7.6ip . (|7.56p . and Tr Ef^ = / 7 , we arrive at (|7.5ip . 

Step IV: The kernel of the trace. We now turn to the task of identifying the null-space of the trace 
operator (17.22[) - ()7.23p . For each k G No we denote by Vk the collection of all vector- valued, complex 
coefficient polynomials of degree < k (and agree that Vk = whenever A; is a negative integer). The 
claim we make at this stage is that the null-space of the operator 

W™' a (Q) \i hl Tr [D~<W}\ G B s J8n) (7.62) 

I J \j\=m—l 

is given by 

V m - 2 + V p m > a (n). (7.63) 

The fact that the null-space of the trace operator (|7.22p - (|7,23p is V p m,a {Q) follows readily from this. 

That (I7.63P is included in the null-space of the operator (|7.62j) is obvious. The opposite inclusion 
amounts to showing that if W G W p '"(O) is such that Tr [D 7 W] = for all 7 G No with \ j\ = m — 1, 
then there exists P m -2 G V m -2 with the property that W — P m -2 G V p m ' a (U). To this end, we note 
that the case m = 1 is a consequence of (|5.20p and consider next the case m = 2, i.e., when 

W G W*' a (Q), Tr [VW] = on 90. (7.64) 

Assume that {Wi}i>i is a sequence of smooth (even polynomial) vector- valued functions in O, ap- 
proximating W in W^' a (0). In particular, 

Tr [VW,-] ^ in L p {dQ) as j -> 00. (7.65) 

If in a neighborhood of a point on 90 the domain is given by {X : X n > cp(X')} for some 
Lipschitz function (p, the following chain rule holds for the gradient of the function Wj : B' 3 X' 1— > 
W,-(X', tp(X')), where is a (n — l)-dimensional ball: 

V Wj (X>) = (VyWjiY'MX')))^ +(jLw j (X',Y n ))\ Yn= ^ Vv(X'). (7.66) 

Since the sequence {wj}j>\ is bounded in L p (B') and Vwj — > in L P (B'), it follows that there 
exists a subsequence {ji}i such that u;^ — > const in L p (B') (see Theorem 1.1.12/2 in [37J). Hence, 
Tr W = Po = const on 90. In view of Tr [W — -Po] = and Tr [VW] = 0, we may conclude that 
W — Po G Vp' a (£~l) by Hardy's inequality. The general case follows in an inductive fashion, by reasoning 
as before with D a W with \a\ = m — 2 in place of W. □ 

We now present a short proof of (jl.8p , based on Proposition 17.31 

Proposition 7.6 Assume that 1 < p < 00, s G (0, 1) and m G N. TTien 

{{il Q lD Q V\ dn }\ a \< m -i ■ V G C °°(R n )} is dense in B™~ l+S (dn) (7.67) 

and 

Wf\\B™~ 1+a (dn) ~ X] [l/allB«(an)) ( 7 - 68 ) 

a|<m— 1 

uniformly for f = {/a}| Q |< m -i £ -B™ _1+s (90). ^4s a consequence, (|1.8p holds. 

Proof. That (I7.67P holds is a direct consequence of (17. 2p and Proposition 17.31 Next, this density 
result and (|7.35p yield the left-pointing inequality in (|7.68|) . As for the opposite inequality, let / = 
{/a}|a|<m-i G S™~ 1+s (90) and, with a := 1 - s - 1/p, consider W := f (/) G H/ p m ' a (0). Then 
Lemma 17.21 implies that, for each multi-index a of length < m — 1, the function / a = Tr [D a Z/] 
belongs to i?p(<90), plus a naturally accompanying norm estimate. □ 

We include one more equivalent characterization of the space B p n ' l+S {d^l), in the spirit of work 
in [1], [49], [60] . To state it, recall that {ej}j is the canonical orthonormal basis in M n . 
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Proposition 7.7 Assume that 1 < p < oo, s G (0, 1) and m E N. Then 



{fa}\a\<m-l G B^- 1+B (dSl) 



' / a G 5|(dfi), Va: |a| < m - 1 
and 

(7.69) 

(^j^fc — Vkdj)f a = Vjf a+ek — Vkfa+ej 

Va : |a| < m — 2, V j, G {1, n}. 



Proof. The left-to-right implication is a consequence of (|7.68p and of the fact that (|7.38p holds for some 
U G W™' a (£l) (cf. Proposition I7,3|) . As for the opposite implication, we multiply {f a }a (as in (|7.11|) ). 
with a function tp G C^°(R n ) which can be assumed to satisfy supp ipndfl = supp^n{(X', X n ) : X n > 
(p(X')} for some Lipschitz function <p : R n_1 — > R. Furthermore, in this latter case, the compatibility 
conditions in (I7,69p become equivalent to (|7.26j) . Thus, given / := {f a }\ a \< m -i whose components are 
as in the right-hand side of (|7.69|) , we may proceed as in the proof of Corollary 17.51 and use the fact 
that f a G B°(dn) for each a to conclude that / G B™~ 1+s (dCl). □ 

7.2 The space of Dirichlet data and the main trace theorem 

In this subsection, we study the mapping properties of the assignment Wwi d k lA/dv k \ , for 

I J 0<fc<m-l 

U G W"jT' a ($7). In order to facilitate the subsequent discussion, for each t G N we consider polynomial 
functions P^f, such that 

n 

V P D <* _ v »dP = H P ^)q^ D ^ Va,/? G N£ : \a\ = = £, (7.70) 

\j\=l-lj,k=l jk 

where d/drjk is the tangential derivative given by 

d d d 

We shall also need Sobolev spaces of order one on dfl. Concretely, let Vtan := (X^ u 3^/^ T 3k)i<k<n 
stand for the tangential gradient on the surface dSl and, for 1 < p < oo, introduce the space 

L£(3J2) := {/ : ||/|Ui (an) := \\f\\ Lp ( 9 n) + l|Vt a „/|| ip(9n) < oo}. (7.72) 

Our main trace/extension result then reads as follows. 

Theorem 7.8 Let Q be a bounded Lipschitz domain in R n and for 1 < p < oo, —l/p<a<l — 1/p, 
s := 1 — a — 1/p G (0, 1) and m G N, consider the mapping 

Tr m _x : W^ a (n) \^\ G L p (dU). (7.73) 

v I OV K ) 0<fc<m-l 

Then its null-space is precisely V™' a (Q), and its image can be characterized as follows. 
Given (g ,g 1 ,...,g m - 1 ) G L p (dtt), set 

/(0.....0) == 90 (7-74) 
and, inductively, having defined {/ 7 }| 7 |<£_i for some £ G {1, ...,m — 1}, consider 

f a :=v«gt + iY: E E ^^M^ Va E N£ : |a| = £, (7.75) 
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where P"j k is any fixed family of polynomials satisfying (|7.70p . Then there exists U G W™ ,a {Q) such 
that Tr m _i U = {gk}o<k<m-i if and only if 

f a £Ll(dn) if \a\<m-2, and f a £B s p (dQ) if \a\ = m - 1. (7.76) 

Furthermore, for Tr m _i as in (|7.73p . set 

W™~ 1+S {dn) := the image of the operator Tr m _i in (17TT31) . (7.77) 

Us llw^-i+^an) := X] ( 7 - 78 ) 

\a\<m—l 



if the families g := {gk}o<k<m-i and {f a }\ a \<m-i are related to oneanother as in (I7,74p - (l7.75p . Then 



i/iis space is independent of the particular choice of polynomials P~fu satisfying ()7.70p and the operator 



Tr m _! : % m ' a (ft) — ► H7~ 1+S (a^) (7.79) 

zs well-defined, bounded, and has a right-inverse. That is, there exists a bounded, linear operator 

Ext : W™- 1+s (dQ) — > % m ' a (0) (7.80) 
smc/i £/ia£ Tr m _i o Ext = /, £/te identity. 

Prior to the proof of Theorem 17.81 we will establish a useful approximation result, extending work 
done in [TJ for the case m = 2. 

Lemma 7.9 Lei Q be a bounded Lipschitz domain in R n with outward unit normal v , and fix 1 < p < 
oo, m G N. assume that {f a }\ a \<m-i is a family of functions satisfying 

f a G Lj(0fi), Va G Nq : |a| < m - 1, (7.81) 

and 

Vjf a+ek - v k f a+ej = Va : ja| < m - 2, V j, fc G {1, n}. (7.82) 

Then there exists a sequence of functions F e G Cq^R™), e > ; such that 

jl«l Tr [D a F £ ] — ► f a in L^dSl) as e -► 0, Va G N£ : |a| < m - 1. (7.83) 

Proof. Since both the hypotheses (|7.81|) - (|7.82|) and the conclusion (|7.83|) are stable under multiplica- 
tion by a smooth function with compact support as in (17. lip , there is no loss of generality in assuming 
that $7 = {X : X n > ip(X')} for some Lipschitz function <p : R™" 1 — > R, and that the functions f a 
have compact support. In this setting, 

[|/a|Ui ( an)~ll/a(-,¥'(-))IUiCR»- 1 )' Va G Nq : \a\ < m — 1, (7.84) 
and the compatibility conditions (|T.82|> can be written in the form 



d 



f a (X',ip(X'))\ = f a+e .(X' :( p(X')) + dj<p{X')f a+en (X f ,<p{X')) 
for a.e. X' G R™" 1 , Va G NfJ : |a| < m - 2, 1 < j < n - 1. 



(7.85) 



Next, fix a nonnegative function rj G Cg°(K n which integrates to one and, for each e > 0, set 
7] £ (X) := e- n r/(X/e). Then, for each |a| < m- 1, e > 0, and X = G R n , consider 

m— 1— 1«| 

F «( X ) : = E M [( X "-^-)) fc / Q+fce „(-^(-))] *%(*')■ (7-86) 

fc=0 
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Based on (|7.85p . a straightforward calculation gives that whenever r := m — 1 — |a| > 1 and 1 < j < 
n — 1, 

djFtix) = ^ +ej (x) + ^^((x n - ¥ p(.)r 1 a^(-)/ a+re „(-^(-))) *ve(x') 

+ L((X n - i p(.)yf a+ren (-,<p(-))) *dj( Ve )(X'). (7.87) 

Thus, if r := m — 1 — |a| > 1, after moving the derivative off of r] E in the last term above we arrive at 
the recurrence formula 



' ^(X) + ^ ((X n - V (-)) r 9 J -(/ ft+f . n (- , *>(•)))) * »fe(Jf') if j<n, 



(7.88) 



Hence, if we now define F e (X) := Ff )(X), an inductive argument based on (17.880 shows that, 

for any multi-index a = (a',a n ) € Nq -1 x N of length \a'\ + a n < m — 1, the difference between 
i\ a \D a F £ (X) and i^(X) can be expressed as a finite, constant coefficient linear combination of terms 
of the type 

e-^[{X n - <p(-)) m - 1 -W- a »d j (f s {.M-))j) * (d^UX'), (7.89) 

where 

l<j<n-l, j3, 7 £ NJ" 1 are such that + + 7 = a', and (5 € Nq, |<5| = m - 1. (7.90) 
Consequently, H7.83P will follow easily once we establish that for every a € Nq of length < m — 1, 

^0, V(0) - U; ¥>(•)) m 4( R ™" 1 ) as £ ^ °i ( 7 - 91 ) 

and that, whenever the indices are as in (|7.90p and X n = 99 (X'), the expression in (17.890 converges to 
zero in Lp(R n_1 ) as e — ► 0. As regards (|7.9ip . we begin by noting that 

m-l-|a| 

i^(X'^(X')) = £ H/ (^)-#')) /.^W^^OM*'-^')^'. (7-92) 
fc=0 ' 1 

Thus, clearly, F^(- , </?(•)) - > /«(• , ¥?(')) m L p (M. n ~ 1 ) as e — > 0, and (|7.91|) is proved as soon as we show 
that V x'Wai' j <f('))] ~^ Vx'[/a(" ! ¥>("))] i n -^p(^™ _1 ) as e ^ 0. We remark that this is obviously true 
if \ot\ = m — 1 (cf. (|7.92p ). so we consider the case when \a\ < m — 2. In this situation, we use (|7.92|) 
to compute, for each k £ {1, ...,n — 1}, 



d 



dX, 



Fl{X'MX')) 



d r 



-1 dY, 



f a+en (Y',<p(Y')) Ve (X'-Y')dY' 



+ 



{dk<p{X') - d k ip{Y'))f a+en {Y\ <p{Y')) rj e (X' - Y') dY 1 



+n{x') 



(7.93) 



with \U(X')\ < C£^| 7 |<m-l I/7O >¥>(•))! *Ve(X'). In particular, ||fc||z, p (R™-i) goes to zero as e -> 0. 



Consequently, gf- F^(X' , tp(X')) 



d 
dX k 



in L p (]R n_1 ) as e -» 0, proving (I7T9TD . 



[/ a (X',^(X'))~ 

Let us now consider the expression (|7.89p when a € Nq with |a| < m — 1 is fixed and the conditions 
in ()7.90p hold. A direct estimate shows that, when X n = cp(X'), the L p -norm of this quantity is 
< Cs m ~\ a \ — > as e —* 0. To finish the proof, take X n = tp(X') in (|7.89|) and, for an arbitrary 
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j G {1, ...,n — 1}, apply d/dXj. Much as we just did, the L p -norm of this quantity is < Ce m_ ' a ' _1 . 
Thus, provided that \a\ < m — 2, this converges to zero as e — > 0. The most delicate case is when 
\a\ = m — 1. In this situation, we write out the terms obtained as a result of making X n = <p(X') in 
(|7.89p and then applying d/dX}, for some fixed k G {1, ...,n — 1}. They are 

i (<p(X f ) - V (Y')yg(Y')(d^r,) £ (X' - Y') dY', (7.94) 



and 

re 



l - r f {^{X') - <p(Y')Y 1 d k (p(X')g(Y')(d 1 ri) £ (X' - Y 1 ) dY' , (7.95) 

where we have set r := m — 1 — \(3\ — a n and g := dj(f$(- , ^(-)))- Above, we have used the fact 
that \a\ = m — 1 forces I7I = r — 1. Our goal is to prove that the L p -norm of the sum between 
(|7.95|) and (|7.94|) . viewed as functions in X E R n , converges to zero as e — > 0. To this end, write 
<p(X') - f{Y') = A(X', Y') \X' - Y'\ + V<p(X') ■ (X' - Y'), where 



i(rt ^.fib»!il). ( , 96) 



then expand 



fax')-v{Y>)) T = ^^(X\Y'f\X'-Yr{V^X')-{X'-Y')) b (7.97) 



a+b=r 



Yl —A(X',Y') a \X'-Y'\ a (Vv(X')) a {X'-Y 



t\(7 



a\a\ 

a+b=r W\=b 



Plugging this back into ()7.94p finally yields 

E E T^i^^y I Mx',Yrg(Y')(e;teM x '- Y ') dY '> ( 7 - 98 ) 

Ci. (J. /TtPn — 1 



a+6=r \"\=b 



where we have used the notation 

Q^ a (X') := (X'Y \X'\ a (d T ri){X') 7 X' € E" _1 , a, r G NJ _1 , a G N . (7.99) 

Each integral above is pointwise dominated by C(||V 1 (p\\L 00 )M.g{X l ) uniformly with respect to e > 
(recall that A4 is the Hardy-Littlewood maximal operator), and converges to zero as e — > whenever 
a > and X' is a differentiability point for the function (p. Thus, since (p is almost everywhere 
differentiable, by a well-known theorem of H. Rademacher, and since M. is bounded on L p if 1 < p < 00, 
Lebesgue's Dominated Convergence Theorem gives that all integrals in (|7.98|) corresponding to a > 
converge to zero in L p (W l ~ 1 ) as e — > 0. 

On the other hand, in the context of (|7.98p . a = forces |<r| = r = (7 + e^j. Note that in general, 
if a = and \a\ = |t|, definition (|7.99p and repeated integrations by parts yield 

f e°>°(x')dx'= [ {x'y{d T r,){x')dx' = {-\)\ T \<y\5 aT , (7.100) 

where 5 aT is the Kronecker symbol. Consequently, as e — > 0, the portion of (I7.98P corresponding to 
a = (and, hence, the entire expression in (|7.98p ) converges in L p (W n ~ 1 ) to 

(-l) r r\(Vip)"< +ek g. (7.101) 
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The analysis of (|7,95p closely parallels that of (|7.94[) . In fact, given the close analogy between (|7.95p 
and (|7.94|) . in order to compute the limit of the former in L p as e —* 0, we only need to make the 
following changes in (|7.10ip : replace 7 + e k by 7, r by r — 1 and then multiply the result by rd^^p. 
The resulting expression is precisely the opposite of (|7.1Uip . and this finishes the proof. □ 

After this preamble, we are now ready to present the 

Proof of Theorem 17.81 The fact that V p m ' a ({l) is the null-space of Tr m _i follows from Proposition ^, 31 
once we notice that (|7.73|) is the composition between (|7.22p and 



-rtm—l+s 



(dn) 3f = {/ a }| Q |< m _! ^ { E ^ va /«} 

\a\=k 



0<fc<m-l 



G L p (dQ), 



(7.102) 



and that the assignment (|7.102p is one-to-one. The latter claim can be justified with the help of the 
identity 

Da = >- ,a, » a ^+ E E^w^"", <«<*> 

\/3\ = \a\ — lj,k=l J 

where p^j! are polynomial functions. Indeed, let / G B p n ^ 1+S (dn) be mapped to zero by the assignment 
([7TT02D and consider U := £{f) G W™' a (VL). Then f a = Tr [D a U\ on dn for each a with \a\ <m-l 
and, granted the current hypotheses, d^/d^ = for k = 0, 1, ...,m — 1. Consequently, (|7.103p and 
induction on \a\ yield that Tr [D a U] = on dn whenever \a\ < m — 1. Thus, ultimately, f a = for 
each a with |a| < m — 1, as desired. In turn, the identity (j7.103H can be proved by writing 



•H D a 



n „ 



■\dxj 

j=l 3 
n otj 



71 n 

n[E&(«. 

i=i k=i 



_d_ 

dxj 



3 <9xi 



II [E /!(«/-/)! & {^d^ ~ dx" 

7=1 i=0 v J 7 k=l 3 



k=l 



CKj-l , 5' 



€=1/ 



(7.104) 



7=1 1=0 y 3 ' k=l K 



-1 x &_ 



5=^ 



and noticing that n?=i Vj 3 d 01 * /di 



,a Q\a\ j Q v \a\ ^ w h ereas (£ k d/dXj — ^jd/dxk) |f= 



-d/<9r. 



Parenthetically, let us point out here that the identity (|7.104p readily proves the existence of some 
polynomial function P"j k such that (|7.70p holds. 

Turning to the characterization of the image of the operator (|7.73p . assume that g k € L p (dQ), 
< k < m — 1, are such that the functions f a defined as in (17.740 - (17. 751) belong to Bp(dQ). The claim 
that we make is that / := {f a }\ a \< m -i G i?™~ 1+s (<9$7) and 



Jfc! 



9k 



E^T* 7 *^' 0<fc<m-l. 



(7.105) 



a|=fc 



Regarding the first part of the claim, by (|7.76p and Proposition 17.71 it suffices to show that (|7.82[) 
holds. We shall prove by induction on i := \a\ G {0, ...,m — 2}. Based on (I7,74p - (l7.75p . we compute 

.0 d 9o 



fa 



" a 9l + E v 



Va : \a\ = 1, 



(7.106) 



1,91=1 



from which the version of (|7.82p with a = (0, 0) is immediate. To prove the induction step, assume 
that (ELE2J) holds whenever |a| < £ - 1. By Lemma E3 there exists F £ G C yo (W l ), e > 0, such that 



\l\<i 



Tr [D^F e ] -> / 7 in Li(dfi) as e -> 0. 



(7.107) 
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From (|7.107p and ([7770]) it follows that for each aeNJJ with \a\ = £ + 1, 

/„ := ^ + Jim £ ^±221 ^ (> Tr [I^PJ - j/* Tr [D?F e ]) in L p (Sn). (7.108) 

\(3\=£+l 

Next, fix an arbitrary a 6 Njj with \a\ = £, choose j, k G {l,...,ra}, and consider the identity (|7. 108|) 
written twice, with a + e& and a + ej, respectively, in place of a. If we multiply the first such identity 
by uj, the second one by Vk and then subtract them from one another, we arrive at 

|/3|=£+1 IX ^ 

By (]7.107p . the above limit is i~ e df a /drjk and this finishes the proof of the induction step. Thus 
(17782]) holds and, as a result, / := {f a }\ a \<m-i G B™- 1+S (dtt), as desired. As for (|7.1U5|) . if we set 

U :=£f £W™' a (Q), (7.110) 
it follows from (17725]) . (17775]) and (177701 that 

f a = v a 9k +J2gi u ^ " Va : |a| = fc, (7.111) 

from which we deduce that i^g/c = v a Y2\/3\=k ffl.^ffi ^ or ea ch multi-index a of length /c. Multiplying 
both sides of this equality by and summing over all a G Ng with \a\ = k finally yields (I7,105p . 

Going further, from (|7.110|) . (|7.105|) and (jl.5j) . we may conclude that {gk}o<k<m-i = Tr m _i^/, 
which proves that the family {gk}o<k<m-i belongs to the image of the mapping f|7.73[> . Conversely, 
if {9k}o<k<m-i = Tr m _i^/ for some function U G W p m,a {Vt), it follows from (|7.70[) and (|7.75p that 
f a = i^Tr[D a U] for \a\ < m - 1. Consequently, f a G B*(dSl) if |a| < m - 1 and f a G L l p {dti) for 
|a| < m — 2, thanks to (|7.69p . This finishes the proof of the fact that (|7.76p characterizes the image of 
the operator (]7.73p . That the space (|7.77|> is independent of the choice of polynomials P"^, satisfying 
(|7.70p is implicit in the above reasoning. Finally, the results in §7.1 imply that the operator (|7.79p is 
bounded. Since as a byproduct of the above proof, the assignment 

B™~ 1+S (dn) b f = {f a }\a\< m -i -> { £ §^/al i e (7.112) 

|a|=fc 

is an isomorphism, we may take Ext in (|7.80p to be the composition between the operator (I7.24p and 
the inverse of the mapping (|7.112|) . This finishes the proof of the theorem. □ 
A specific implementation of the algorithm (I7.74p - (l7.75p is discussed below. 

Corollary 7.10 Assume that Q be a bounded Lipschitz domain in M n and fix 1 < p < oo and — 1/p < 
a < l-l/p, s := 1 — a — 1/p G (0,1), m G N. For a family (g , g 1 , g m -i) G L p (dtt) set /( 0) ...,o) '■= 9o 
and, inductively, if {/ 7 }| 7 |<^_i Ziaue already been defined for some £ G {1, ...,m — 1}, sei 

fa ■= v a 9t + f £ ^ ( V t- W;, V a G Nq : |a| = £, (7.113) 



W=\s\ 



Hi 



where (-)j is the j-th component. Then g = (go,gi, ...,g m — i) belongs to W™ +s {dVL) if and only if 



f '■= {fa}\a\<m-i belongs to B™ 1+s (dQ), in which case \7.105\) also holds. 
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Proof. For any two multi-indices a, (3 G Nq of length t written as a = + • • • ej e and (3 = e kl + ■ • • e kt , 
a direct calculation yields 

a 

vPn* _ ^ = . . . v kl _ r _ lVk _ r+1 ■ ■ ■ v k - d n ■ ■ ■ d jM d k£ _ r+1 ■ ■ ■ d H . (7.114) 

r=0 UT k e _ r j e _ r 

In order to be able to re-write (|7.114p in multi-index notation, it is convenient to symmetrize the 
right-hand side of this identity by adding up all its versions obtained by permuting the indices j\ , ji 
and ki,...,kt. In this fashion, we obtain 

^-^-i^V V (*-r-l)l rl (/-r-l)l rl + , d B 

r=0 n+6+ej=at, \S\=r n 1 

1+e+e k =f3, \B\=r 

This is a particular version of (|7.7U|) . where the intervening polynomials are identified explicitly. If we 
now implement the algorithm (I7.74| )-( j7.75p . for each a G Nq with \a\ = £ we arrive at 

f*- v nt + Tx l^^ 1, — ~ { — ji — ~\ — oi v (7 - 116) 

\{3\=£r=0 n+6+ej=ct,\S\=r ^ ' K J 

■y+0+e k =0, \$\=r 

Next, we replace (3 by 7 + 9 + e&, eliminating the sum over and make use of the identities 



(*_-r-l)!.j*_ n df 
in order to transform (I7.116P into 



£ {1 ^P h ^ = ^ E^&.=^fh (7-117) 



a! ^— i, v-^ (I — r — 1)! r! r! 



= + £ ^ (7-118) 



\B\ = \S\=r 



which is equivalent to (|7.113p . □ 

The space (|7.77p takes a particularly simple form when m = 2. Indeed, as a direct consequence of 
(|7,113p in which we take t = 1 we have: 

Corollary 7.11 For each Lipschitz domain £1 C M. n and each 1 < p < 00 , s G (0, 1) , 

wJ+'Can) = {(50,51) € Ll(dQ) © L p (an) : ^ + v tan 5o e £pW>- (7.119) 

This has been conjectured to hold (when s = 1 — 1 /p) by A. Buffa and G. Geymonat on p. 703 of [8 J . 

Finally, we comment on how (|7.77p relates to more classical spaces of higher order traces when 
Q C R n has a smoother boundary than mere Lipschitz. Specifically, fix m G N, m > 2, p G (l,oo), 
s G (0,1) and assume that dQ is locally given by graphs of Lipschitz function <p : R™" 1 — > R with the 
additional property that V</? belongs to M£>™ _2+s (R n_1 ), the space of (pointwise) multipliers for the 
Besov space i?™ _2+s (R™ _1 ) (cf. [39], [41j ). Then, for each non-integer 1 < /x < to — 1 + s, one can 
coherently define the space Bp(dft) by starting from (R n_1 ) and then transporting it to dfl via a 
smooth partition of unity argument and by locally flattening the boundary. In fact, we arrive at the 
same space by taking the image of the trace operator on dfl, acting from i?p +1 ^ p (R n ). 
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Proposition 7.12 Assume that p G (l,oo), s G (0,1) and that C M n is a Lipschitz domain whose 
boundary is locally described by means of graphs of real-valued functions in M n_1 whose gradients 
belong to M5 m - 2+s (M n ^). Then 



m—l 

w™~ 1+s (dn) = Yl B™- l - k+s (dn). (7.120) 

fc=0 

In particular, this is the case if dtl G C m ~ 1 ' for some 9 > s. 

Proof. In one direction, (17.690 and lifting theorems imply that if {fa}\a\<m—i 

G B™~ 1+S (dn) then 

f a G S^~ 1-|Q!|+S (90) for each a with \a\ < m - 1. Hence, g k := £ H=fc § u a f a € B™- 1- ^^) for 
each G {0, ...,m — 1}, so the left-to-right inclusion in (|7.120p follows from the fact that (|7.112p is an 
isomorphism. 

As for the opposite implication, given {gk}o<k<m-i £ ©fc=o B™~ 1 ~ k+S (dtt) , define {/a}|al<m-l as 
in (|7.74|) - (|7.75|) . Granted the current assumptions on 50, an argument based on induction and the fact 
that v G MB r p n ~ 2+s (dn) shows that f a G S^" 1_|a|+S (50) for each |a| < m - 1. In particular, (|T7?6j) 
holds which proves that {gk}o<k<m-i £ W™ +B (dtl). This shows that the right-to-left inclusion in 
(|7.120j) also holds, thus finishing the proof of the proposition. □ 

8 Proof of the main result 

8.1 The inhomogeneous Dirichlet problem 

Theorem 11.11 is a particular case of Theorem 18. 1| concerning the solvability of the inhomogeneous 
Dirichlet problem 

' A(X,D X )U = T inO, 

d k lA (^'1) 
—r = g k on dQ, < k < m - 1, 

in the space W p ,a (J2). Note that for any operator A as in §6.1 we have 

A(X, D x ) : W™' a {Q) — » F- W ' a (fi) (8.2) 

boundedly. Thus, granted the membership of U solving (|8.ip to W™' a (0), it follows from Theorem 17.81 
that necessarily T G V p ' m,a (Q) and 5 := {gfc}o<fc<m-i € W"™ -1+s (<90). Moreover, 

lis II w-™-^ ( 5 n) + ll-^lly-^tn) - c 'll w llw^ l ' a (n)- ( 8 - 3 ) 
The converse direction makes the object of the theorem below. 

Theorem 8.1 Let Q be a bounded Lipschitz domain in W 1 and assume that the operator A is as in 
§6.1. Then there exists c > such that if (|1.13p is satisfied then the Dirichlet problem (|8.ip has a 
unique solution U G W™' a (Q) for any given T G V p m ' a (Q,) and g := {gk}o<k<m-i £ W™ +s (dQ). 
Furthermore, there exists C = C(dQ,A,p,s) > such that 

\M\w™- a (n) < c (\\ a \\w; n - 1+s (dn) + H^lly-" 1 --^))- ( 8 - 4 ) 

Proof. We seek a solution for (|8.ip in the form U = Ext (g) + W, where Ext denotes the exten- 
sion operator from Theorem 17.81 and W G V p m,a (£l). Note that, by Theorem 17.81 this membership 
automatically entails d k W jdv k = on dQ for k = 0, 1, ...,m — 1, so it suffices to take 

W := A(X, Dx)' 1 (T - A(X, D x )Ext (g)) G V p m ' a (fl) (8.5) 
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which, by (|8.2p and Theorem 16,31 is meaningful. As for uniqueness, let U G W p ' (Q) solve (|8.ip with 
J 7 = and = 0, < k < m — 1. Then the function W belongs to V^™' ^), thanks to Theorem 17,81 
and is a null-solution of A(X,Dx)- In turn, Theorem 16.31 gives that U = 0, as desired. Finally, (18, 4p 
is a consequence of the results in §7. □ 
We conclude this subsection with a remark pertaining to the presence of lower order terms. More 
specifically, granted Theorem 18 .1\ a standard perturbation argument (cf., e.g., |27j ) proves the follow- 
ing. Assume that 

A{X,D X )U:= Yl D a {A aP (X)D l3 U), X G 0, (8.6) 
0<H,|/3|<m 

where the top part of A{X,Dx) satisfies the hypotheses made in Theorem 11.11 and the lower order 
terms are bounded. Then, assuming that either (18. 9p or (jl . 131) holds, the Dirichlet problem (18. ip is 
Fredholm with index zero, in the sense that the operator 

W p m ' a (tt) (A(X,D X )U , {d k U/du k } < k < m ^ € V- m ' a (Q) W™- 1+S (dtt) (8.7) 

is so. Furthermore, the estimate 

\\U\\w™' a (n) < c (WfWv-^^n) + WdWw^-'+^dn) + \\ u h p (n)) (8.8) 
holds for any solution U 6 W™' a {Q) of IJOl . 

8.2 Further comments and the sharpness of Theorem 18.11 

A byproduct of our proof of Theorem 18.11 is the following. Assume that C W 1 , jC are as in the first 
paragraph of the statement of Theorem 18. 11 Then there exists e > 0, depending only on the Loo-norm 
of the coefficients and the ellipticity constant of C, with the property that the Dirichlet problem (|8.ip 
with data from W™~ l+s {dQ) has a unique solution in W™' a (fl) granted that 

|2 -1 -p -1 | <e and |o| < e. (8.9) 

To justify this claim, we rely on Theorem 17.81 and, using the same strategy as before, reduce matters 
to proving that the operator (I6.28P is an isomorphism. When a = and p = 2, our assumptions on 
A(X, Dx) and the classical Lax-Milgram lemma ensure that this is indeed the case. Then the stability 
theory from |32j, |52j allows us to perturb this result, i.e., conclude that (16.281) is an isomorphism 
whenever (|8.9p holds, as soon as we show that the scale V^ m ' a (f2) is stable under complex interpolation. 
That is, if 1 < pi < oo, -1/pi < en < 1 - l/pi, i G {0, 1}, 9 G (0, 1), 1/p = (1 - 9)/p + 9/pi and 
a = (1 — 9)clq + 9ai, then 

[v^ ao (n),v p ^(n)] e = v™> a (n), (8.10) 

where [•, -]q denotes the usual complex interpolation bracket. In the proof of (|8.10p we may assume 
that £1 is a special Lipschitz domain and, further, that £1 = , by making the change of variables 
described in §5.2-§5.3. In this latter setting, it will be useful to note that 

[L po (Rl, x a n 0P0 dx),L pi (Rl, x a n lPl dx)] 6 = L p (Rl, x * dx), (8.11) 

granted that the indices involved are as before, which follows from well-known interpolation results for 
Lebesgue spaces with change of measure (cf. Theorem 5.5.3 on p. 120 in [6]). Then (I8.10|) follows easily 
from (|8.1ip . the fact that for each a G Nq with \a\ = m the operator D a maps the scale V p ' a (Wt ) 
boundedly into the scale L p (M+, x^ dx), and (|4.23|) - (j4.25j) . This finishes the proof of the claim made 
at the beginning of this subsection. 
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In turn, the aforementioned result can be viewed as an extension of a well-known theorem of 
N.Meyers, who has treated the case m = 1, I = 1 in |4*2|. The example given in §5 of |42j shows 
that the membership of p to a small neighborhood of 2 is a necessary condition, even when d£l is 
smooth, if the coefficients A a p are merely bounded. For higher order operators we make use of an 
example originally due to V.G.Maz'ya |36j (cf. also the contemporary article by E.De Giorgi |16j). 
Specifically, when m G N is even, consider the divergence-form equation 

Ai m-1 £ 4 A^^U = in !!:={IeK": \X\ < 1}, (8.12) 

where £4 is the fourth order operator 

C 4 (X,D X )U := a^ 2 U + bj2^{T^^ d j U )+ h iZ did 3^T^^ U 

i,j=l i,j=l 

+c £ ^( ^ ^a). (8.13) 

i ,j,k,l=l 



Obviously, the coefficients of C±(X, Dx) are bounded, and if the parameters a, b, c G M, a > 0, are 
chosen such that 6 2 < ac then C along with A2 m_1 £ 4 A2" 1 "" 1 are strongly elliptic. Now, if VKp 
denotes the usual L p -based Sobolev space of order s, it has been observed in |36] that the function 
U{X) := \X\ e+m ~ 2 G W 2 m (0) has TrW G C°°(<90) and is a weak solution of (I8TT2D for the choice 



g:=2-^W^- (n - 1)( f n + c) . (8.14) 
2V4 a + 26 + c v; 

Thus, if a := (n — 2) 2 + e, 6 := n(n — 2), c := n 2 , e > 0, the strong ellipticity condition is satisfied 
and 9 = 9(e) becomes 2 - n/2 + ne 1 / 2 /2y/4(n - l) 2 + e. However, U G WZ"(JT) if and only if 
p < n/(2 — 0(e)), and the bound n/(2 — 9(e)) approaches 2 when e — > 0. An analogous example can 
be produced when m > 1 is odd, starting with a sixth order operator £q(X, Dx) from |36j. In the 
above context, given that W„(0) VMO(f2), it is significant to point out that both for the example 
in |42j, when n = 2, and for (|8.12[) when n > 3, the coefficients have their gradients in weak-L n yet 
they fail to belong to W^(f2). 

Of course, condition (|1.3|) ensures that the left-hand side of (|1.13|) is always finite but it is its 
actual size which determines whether for a given pair of indices s, p, the problem (|1.2p . (|1.6I) . (|1.9f) is 
well-posed. Note that the maximum value that the right-hand side of (|1.13j) takes for < s < 1 and 
1 < p < 00 occurs precisely when p = 2 and a := 1 — s — 1/p = 0. As (]l,13p shows, the set of pairs 
(s, 1/p) G (0, 1) x (0, 1) for which (|8.ip is well-posed in the context of Theorem 18. II exhausts the entire 
square (0, 1) x (0, 1) as the distance from v and the A a ^s to VMO tends to zero (while the Lipschitz 
constant of O and the ellipticity constant of C stay bounded). That the geometry of the Lipschitz 
domain f2 intervenes in this process through a condition such as (|1.13p confirms a conjecture made by 
P. Auscher and M. Qafsaoui in [5] . 

While the main aim of the present work is the consideration of higher-order operators with coef- 
ficients in Lqo, Theorem 18.11 (and, with it, Theorem II .ip is new even in the case when m = 1 and 
A a g G C lxl (i.e., for second order, constant coefficient systems). It provides a complete answer to the 
issue of well-posedness of the problem (I8.ip in the sense that the small mean oscillation condition, 
depending on p and s, is in the nature of best possible if one insists on allowing arbitrary indices p 
and s. This can be seen by considering the following Dirichlet problem for the Laplacian in a domain 
n C R n : 

AU = in n, TrU = ge B'(d£l), \U\ + \VU\ G L P (Q, p(X) p{1 - s) - 1 dX). (8.15) 

It has long been known that, already in the case when d£l exhibits one cone- like singularity, the 
well-posedness of (|8.15p prevents (s, 1/p) from being an arbitrary point in (0, 1) x (0, 1). At a more 
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sophisticated level, the work of D. Jerison and C.Kenig in the 1990's shows that (|8.15p is well-posed 
in an arbitrary, given Lipschitz domain f2 if and only if the point (s, 1/p) belongs to a certain open 
subregion of (0, 1) x (0, 1), determined exclusively by the geometry of the domain £2 (cf. |29j). 
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